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TRANSIENT AND STEADY-STATE PROCESSES IN AN AUTOMATIC 
RANGE INDICATOR 


PART I. DESCRIPTION OF THE DEVICE'S FUNCTIONING AND USE OF THE APPARATUS OF 
STEP FUNCTIONS IN THE ANALYSIS OF INDICATOR PROCESSES 


F. M. Kilin 
(Leningrad) 


The transient and steady-state processes in an automatic range 
indicator are considered, account being taken of the specific idio- 
syncrasies characteristic of its functioning, Among these are the dis- 
continuous character of the processes in the individual elements of the 
indicator, the variability of the design parameters which change by 
jumps in correspondence with the applied pulses, and the temporal 
transformations of pulses in the coincidence amplifiers, 


1, Introduction 


Automatic range indicators, conventionally called auto-range-finders (avtodal'nomer), are among the compo- 
nents of modern radar and radio-navigation systems. In these systems the auto-range-finder may execute various 
functions, We shall consider these functions in terms of an example, the auto-range-finder of a radar set, 


Figure 1,a shows what we shall call reference 
pulses, formed by circuits of the radar transmitter. 


Figure 1,b shows pulses taken off the output of 
this radar's receiver, It follows from this display that 
there are some objects in the space scanned by this 
radar. In this case, an auto-range-finder can be used 
to solve the following problems, 


1. Automatic tracking by the return pulses of 
any one object, for example, tracking on the A pulses 
(Fig. 1,b). 


2. Automatic selection of signal pulses, If the 
auto-range-finder is tracking on the A pulses then, 
by means of this, it is possible to isolate the A pulses 
Fig. 1. and to apply them to the circuits following the auto~- 


range~finder (Fig. 1,c). 
3. Automatic measurement of the time interval between the reference pulses and one of the return pulses 


on which the set is tracking. In the case of the auto-range-finder considered, the time interval measured would 
be that denoted by t on Fig. 1. 
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2. Auto-Range-Finder Block Schematic 


. Figure 2 gives a simplified block schematic of a radar auto~range-finder. The auto-range-finder consists 
of two parts, a tracking part and a selection part. Included in the tracking part are a time discriminator TD, an 
operational amplifier OA, a regulated delay circuit RDC and a tracking pulse generator TPG. 


.-Selection portion of auto-range finder 


R of +} 
RDC 
T ai, 
2 
CA-1] (CA-2 
DD 
OA 
- ol 
Output of the 
| CDC SPG CA selection 
Output of the 


| Selection portion of auto-range finder selection 


Fig. 2. Block schematic of an auto-range-finder device: A is the antenna, 
AS is the antenna switcher, R is the receiver, T is the transmitter and SD is 


the scope display. 


Figure 3 provides graphs delineating the functioning of the auto-range-finder. Graphs 3,a represents the 
sequence of reference pulses with period T). These impulses are applied at the input of the RDC from the trans- 
mitter, The action of the reference pulses on the RDC causes the latter to form its own pulses, shown in Graph 3,b, 
These pulses are lagged in time by a certain amount with respect to the reference pulses. The magnitude of the 
pulse lag is proportional to the signal applied to the RDC from the operational amplifier. From the RDC, the 
lagged pulses are applied to the input of stage I of the TPG, and trigger it. As a result of this, stage I of the TPG 
forms pulses of height Ey and length oy, the leading edges of which coincide with the leading edges of the cor- 
responding lagged pulses, In stage II of the TPG are formed pulses of the same height, EO, and length, ap», as the 
pulses formed by stage I, but displaced from the latter by the quantity a». Triggering of stage II of the TPG is 
effected by the trailing edges of the stage I pulses. 


The sequences of pulses formed by stages I and Il of the TPG are described by the expressions, respectively, 
Ey = Eo Y(t), Eqy = Ep Ye(t). The graphs of y,(t) and y,(t) are given in Fig. 3. Graph 3,f shows the signal pulses 
from the receiver output which, in what follows, will be denoted by U,. A short description of the individual 
elements in the auto-range-finder is given below. 


The time discriminator generates a voltage (or current) proportional to the temporal displacement of the 
pulses of one sequence with respect to the pulses of another sequence, In the auto-range-finder design, the TD 
is a sensitive element by means of which may be detected the temporal shift of the pulses reflected from the 
target with respect to the tracking pulses, The principle of operation of the TD may be explained by the simpli- 
fied circuit shown in Fig. 2. This circuit consists of two coincidence amplifiers Ca-1 and CA-2 and a differential 
detector. At the input terminals of the TD (1 and 1") are applied pulses U,(t) from the receiver, which we shall 
call return pulses, From the output terminals 4 and 4’, is taken a voltage proportional to the temporal displace- 
ment of the return pulses with respect to the tracking pulses, The TD has additional terminals 2 and 2", 3 and 3°, 
which receive pulses from the tracking pulse generator corresponding to Ej and E;,. The functioning of the TD 
consists of the following. 


Each coincidence amplifier of the TD passes those portions of the target return pulses which coincide, in time, 
with the corresponding pulses applied to this coincidence amplifier from the TPG. Consequently, the cross-hatched 


ly, 


portion of pulse A, of length ay (Fig. 3,f), is applied to coincidence amplifier CA~1, and the remaining portion 


of this pulse, of length ag, is applied to coincidence amplifier CA-2, Both coincidence amplifiers have identical 
circuits and parameters, 


Figure 3,g shows the signal pulses, U,(t), after their temporal transformation in the coincidence amplifiers, 
The gain of the coincidence amplifiers is denoted by Ky. The quantity Aa[n] is the temporal shift of the n'th 
target pulse with respect to the center line of the tracking pulses 
with the same ordinal number. The pulse transformations are 
described by the expression KylJ,(t)y(t), where y(t) = 74(t) 
(Fig. 3,e). Thus, in the passage of the return pulses from the 
= % ; target through the coincidence amplifiers, two sequences of 
7 RDC lagged pulses pulses are formed and then applied to the differential detector 

DD (Fig. 2). In the differential detector, the difference of the 
¥(t) constant components of these two pulse sequences is isolated. This 

of difference is proportional to the magnitude of Aafn]. 


| AC = The operational amplifier amplifies the voltage formed in 
S the circuits of the time discriminator, and transforms this voltage 

i cee in-accordance with some law described by the corresponding dif- 
€ ferential equation, The output voltage of the OA determines the 
magnitude of the pulse delay in the RDC, 


The functions executed by the RDC and the TPG have been 
discussed above. Included in the selection portion of the auto- . 
range-finder are a constant-delay circuit, CDC, a selection pulse 
Keo) oa generatar, SPG, and a coincidence amplifier, CA (Fig. 2). The 

selection portion of the design is controlled by the delayed pulses 
oe of the CDC. Pulses are formed in the SPG which coincide in 
8 [n} iF time with the pulses from that one of the targets which the auto- 
range-finder is supposed to track. The selector pulses, in turn, 
— control the functioning of the coincidence amplifier, CA, Asa 
result, of all the pulses applied to the CA from the receiver out- 
Pig. 3. put, only those pulses are passed which are returns from the one 
target being tracked, 


Return pulses 


The temporal placement of a selector pulse with respect to the corresponding transmitter reference pulse 
depends on the magnitude of the signal at the OA output, i.e., it is continuously corrected by the tracking portion 
of the auto-range-finder. Therefore, the selection portion of the auto-range-finder, working in conjunction with 
the tracking portion, is conventionally called an automatic selector. 
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Descriptions of various designs of time discriminators, differential detectors and coincidence amplifiers are 
to be found in the literature [1-4]. For the regulated delay circuit in the auto-range-finder, wide use is made of 
phantastron circuits, phase-inverters of various types and other devices, Works [1, 3, 5] describe individual auto- 


tange-finder designs, 


3. The Auto-Range Finder as a Servo System 


For analyzing the functioning of an auto-range-finder, it is convenient to estimate the temporal placement 
of the return and tracking pulses with respect to the corresponding reference pulses, Figure 3 showed the reference 
pulses with the ordinal number n, The center of a reference pulse is taken as the origin of coordinates. In this 
coordinate system, the position of the center of a return A pulse is determined by the magnitude of %,[n], 
which varies from one pulse to the next in correspondence with the variation in the ordinal number of the reference 
pulse, 


On the same figure there is shown the position of the median line of the corresponding pair of tracking 
pulses, which is estimated by the magnitude of @{n]. If the auto-range-finder carried out ideal tracking, the 
median line of the tracking pulses would be exactly superposed on the median line of the return A pulses, Actually, 
the functioning of a real circuit is characterized by a tracking error, which is defined by the expression: 


Aa [n] = 9, [n] — 9 [nj. 
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The quantities $[n], and Acf[n] are step 
functions, The independent variable n takes only inte- 


TD sex (nf OA gral values; n= 0,1, 2,... A definition of step functions 
3 is given in [6]. In what follows we shall enclose the 
(al independent variable n of a step function in square 
brackets. 
TPG RDC 
To render more graphic our consideration of an 


auto-range~finder as a servo system, we shall use as a base 
the circuit shown in Fig. 4. According to the circuit 
diagram, two quantities ®[n] and $,[n] are applied to 
the time discriminator (TD). In practice, the processes in the time discriminator, described above, reduce to a 
comparison of the magnitudes of $[n] and $,[n], as a result of which an error signal, proportional to Aaf{n], is 
formed at the output of the TD. The error signal is then applied to the input of the operational amplifier (OA) 

in which it is transformed to the quantity U,. In its turn, the signal U, acts on the RDC, engendering a propor- 
tional variation of the quantity 9. Such control of the quantity $ amounts to this latter repeating all the varia- 
tions of the quantity 9. 


Fig. 4. Block schematic of an auto-range-finder as 
a servo system. 


4. Auto-Range-Finder Equations for Noiseless Input to the Time Discriminator 


In this paper we consider a working regimen of the auto-range-finder in which the return pulses from the 
target being tracked coincide in time both with the first and the second tracking pulse (Fig. 3), With such an 
orientation of the return pulses, these latter are divided into two parts by the time discriminator. For investi- 
gating auto-range-finder dynamics, we use the methodology presented in [7] which allows the piece-wise- 
constant variability of the parameters in the individual elements of the device to be taken into account. This 
methodology is a further development of the apparatus of step functions, the fundamentals of which are presented 


in [6]. 


We mention in passing, that the apparatus for investigating the processes in pulse systems of automatic 
control is presented in [8]. However, the methodology described in this work is only applicable to systems with 
constant parameters, Moreover, the questions broached in this work are treated less completely than in [6]. 


We consider the functioning of an auto-range~-finder in the case when the pulses reflected from the target 
are rectangular, or nearly so, With this assumption, the sequence of pulses applied to the input of the time dis- 
criminator can be described by the expressions: 


Us; (t) (¢), 
= [H (t —kT) — H (t —a—AT)}. 
k=0 


Here, H(t) is a function defined by the equations: 


0 for t<1, 


¥ oft) is a pulse function of unit height E., T and a are constant magnitudes of, respectively, the height, repeti- 
tion period and length of the pulses reflected from the target (Fig. 5). 


The tracking pulses, formed by stages I and II of the TPG, are described by pulse functions of the form: 


= + — KT) — H(t —A% — 


= — — AT) — H(t — Ad — 


k=0 


49, is the magnitude of the variable pulse lag in the RDC. The pulse functions U,(t) and y(t) = y4(t) —Ye(t), 
are shown on Fig. 5. 


n net 


| 


Fig. 5. Transformation of the signal pulses in the time discri- 


minator, 


Time discriminator equation, The coincidence amplifiers in the time discriminator implement the function 
of high-speed switching. In practice, their transient responses do not affect the magnitude of the output signals. 
Thus, the processes in the time discriminator reduce to the processes in the circuit shown in Fig. 6, In this circuit, 
switch Ky carries out the function of coincidence amplifier 
CA-1, closing periodically for a length of time equal to the 
length of tracking pulses 1, Switch Kg carries out the function 
of coincidence amplifier CA~-2, also closing periodically but 


! 
4 


CA) 


ly (7D , during intervals of time equal in length to tracking pulses Il. 
bo =. The voltages from the terminals of generators E; and E, are 
T | determined by the formula: 
t E, = —E, = KoE y(t), 
P 7 { where Ky is the gain of the coincidence amplifiers. 


Differential vector Thus, both generators supply voltages which are equal 
2 in absolute value but opposite in polarity, The functioning 
_ ee _— of the time discriminator, in accordance with the circuit 
. of Fig. 6, reduces to the following processes. 


With switch Ky closed, a portion of pulse U, with the ordinal number n (Fig. 5), of length a), and coinci~ 
ding with tracking pulses I, is passed by, and charges condenser Cy. The remaining portion of this pulse, of 
length cgn and coinciding with tracking pulses II, is passed with negative polarity when switch Kg is closed, With 
this. condenser Cy is discharged. As a result of the action of both pulses on condenser CO, a voltage is formed, 
Proportional to the difference of the pulse lengths a}, —cgn. The lengths a,,, and cgn, are related to the auto~ 
range~-finder error in the following manner: 


(1) 
Zin — = 2Aa [nm]. 
Based on what has been presented, the equation of the time discriminator is written as: 
[1 (t) D + 1] (t) = Kok ho (t) 1 (2). (2) 


Here, D = d/dt is the differentiation operator and 1, is a time constant which varies in accordance with the 
expression: 


(t) = (t) + (4), = A(t) — (t). 


The magnitude of T; varies in correspondence with the pulse function U,(t). This variation arises from the 
circumstance that, if there are input pulses, one of the diodes of the differential detector is conducting and, in 
correspondence with this, the circuit has one set of parameters, while with no input pulses, the diodes become non- 
conducting and the circuit assumes another set of parameters. 


In accordance with the circuit of Fig. 6, the time constants Ty and Tyy are defined by the parameters of 
the differential detector in the following way: 


= = CoRo, R= + fi, 
where R, is the internal impedance of diodes Ly and L, (Fig. 6). If Rp >> Ry, the formula for Ty assumes a simpler 
form: 
CoR,. 


Operational amplifier equation, In the general case, this equation is written in the following form: 


+ Ux (t) = 0, =2, 3,..., 7) (3) 


k=1 
where 6 ,, is a Kronecker delta, 


Equations of the regulated delay circuit (RDC). In the RDC, pulses are delayed by an amount proportional 
to the signal taken off the output of the operational amplifier. This can be expressed by the following equations: 


= A$, Ad = (4) 
_ Here, 9 is the constant delay of the RDC, A® is the variable lag of the RDC and o is the constant coef- 
ficient of proportionality. 
Auto-range~finder error equation, In paragraph 3 the auto-range-finder error equation was obtained in the 
form: 


Aa[n}] [n] — 
With Equation (4) taken into account, the auto-range~finder error equation can be put in the form: 
Aa [n] + Ad [n] = AS, [n}, A%, [n] = 9; [n] — %. 


5. Transformation of the Equations of the Time Discriminator and the Operational 


Amplifier 


In accordance with (2) and (3), the processes in the time discriminator and the operational amplifier are 
described by the equations: 


[D + ay (t) 


Dd) Ue (t)=0, =2, 3,..., 7). (6) 


k=l 


Tio Ti 


(t) = (t) + (t), an = 


fr (t) = 


n= 


To shorten the writing and the mathematical computations, we shall employ matrix notation, We write 
Equations (6) in matrix form 


[DI + a(t) U (t) = F(t). (7) 
Here, 1 = (J, k = 1, is the unit matrix, 
a(t) = ho (t) + a", 
a? = | ax |, a’ = | (j, 2, 


The functions U(t) and F(t) are column matrices of the form: 


(8) 


U, (t) 0 
vw=|...]. Fol... 
U,(t) 


The symbols 0 and 1 in Expressions (8) are indicative of the fact tliat the coefficients correspond to the para- 
meters of the system in the time intervals determined by the pulse functions p(t) and p(t). The matrix form of 
writing the equation in (7) is more general than the form used in (6), since the case is also included when the para~ 
meters of the operational amplifier have a piece-wise constant form of variability. 


For the matrices, DI + a’ and Di +a", we find the adjunct matrices, denoted by 
A°(D) =| A*(D) = | (D) 
and the determinants: 
A (DI.+'a°)'= [DI + A (DI + a2) =| DI +4341. 


We now determine the function U(t) in the interval nT < ts (n+1)T. This interval is decomposed into 
two regions, shown on Fig. 5. By solving Equation (7), we obtain the following expression for U(t):* 


in region I 
U (t) = tn} + — nT) — h(t — 


Tio (9) 
h(t)=0 for t<0, nT <t<¢nT +a; 
in region II 
KE, 
U (t) = [n} + (a) — 2h (a — (10) 
nT a<t<(n+1)T. 
Here, 
U [n] = {U (t)}:-nr- (11) 
The matrices owt and oat have the following expressions: 


* The solutions of (9) and (10) may be obtained by the methods described in [9] in pp. 99-104, 


In turn, the functions B YK) and Bic!) are found from the operational representation: 


(P) 


Dol +e) (s=0, 1). 


L [Bix (t)] = 


In certain cases it is more convenient to express the matrices e* * and eat in another form, In particular, 


they may be represented by series (Cf. [10], pp. 85-89): 


+Fe_..., 


The function h(t) is a column matrix of the form: 


defined by the expressions: 


t 
h(t) (x) dz, fo(t) = 
0 


or, in expanded form, taking (12) into account: 


t 
hy (t) = | dee 
0 


The elements, hj(t), could also be found directly from the expressions 


(P) 
L [hj (t)] = 


where L is the operator of taking the Laplace transform. 
By setting t = (n + 1)T in Equation (10), and introducing the nomenclature 


b= (T—a)eg—o%e 
= e-(T—2) [h (a) — h(a — 


we obtain the recursion relationship 
U [n +1] [nn] =§[n], 


| 
hy (t) 
h(t)= 
| h;(t) 
| H (t) 
0 
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from which the step function U[n] may be determined, 


The step function U[n] is related to the function U(t) by Equation (11), i.e., it characterizes the variation 
of U(t) at the discrete moments of time, 0, T, 2T,.... In the intervals between these moments, the variation 
of U(t) is determined, for the respective regions, by Expressions (9) and (10), if we substitute in these equations 
the values of the step function U[n] which satisfy Equation (15), In the sequel, we shall call Equation (15) the 
discrete equation of the time discriminator of the operational amplifier. 


The step function U[n] gives a rough representation of the character of the transient and steady-state pro- 
cesses in the system investigated, However, such a representation of the system's dynamic properties is completely 
sufficient for the solution of certain applied problems, For the investigation of processes in a system by means of 
the step function Uf[n], the initial Equation (7), with variable coefficients, is replaced by Equation (15), with 
constant coefficients, the solution and investigation of which pose no particular difficulty, The operational cal- 
culus methods for solving equations of the type of (15) are well developed, and are almost identical with the 
analogous method of solving differential equations with constant coefficients, Thanks to this, the analysis of pulse 
systems with piece-wise constant variability is simplified. 


By means of the time discriminator and the operational amplifier, the retention and transformation into 
continuous signals proportional to &a of the pulse signals are implemented. Such a type of transformation will 
be carried out in the given system if the matrices e~#°a and e~* * are amenable to the approximate representa- 
tions: 
with this, 


| <1. (17) 


Here, the signs | | tho, values of the of toh clement Sy a, 
and of the element atk by the repetition period T. 


By substituting, tn the expression for h(t) in (14), the matrix e-2%t —T) in the form of a series, as in (13), 
and by then carrying out the integration, we obtain: 


If Conditions (17) hold, we have 


h(a)=foe, fo= 


Equations (18) and (1) give us: 
h(a) — 2h (a — an) = 2Aa [n] Jo. 
It follows from this that: 


= Aa [n] 


Thus, Equation (15) can be rewritten in the form: 


2K oli, 


U [n + 14] — bU [n] = oF Aa [n] fo- (19) 


According to (19), the form of the step function U[n] depends on the initial conditions for n = 0 and on the 
form of the step function Aa[n], which expresses the variation of the tracking error from one pulse to the next, 
Thanks to this, signals proportional to Aaf[n] are formed in the elements of the given system. 


6. Transient and Steady-State Processes in the Time Discriminator and Operational 
Amplifier 


If the conditions in (17) hold, the expression for the coefficient b can be simplified, and given by the expres- 


sion: 

(20) 

=e-aT 

aa + a! (T — a) 
a= 
If we substitute the value of b in Equation (19), we get 

U [n + 1] — [n] = Aa [n] fo. (21) 


Tlo 


If we denote by U[0] the initial condition for U[n] for n = 0, we find from Equation (21), the sequence U[1], 
U[2], ..., U[n]. The expression for U[n] is written in the form: 


U [n] = [0] + 


(n—K—1) Aw ho. (22) 


k= 


The matrix e~** has the expression 
= Bix (t)) 4, 2, ..., (23) 
The functions By) are found by the operational transformations: 


(P) 
L [Bix (1 = (24) 


where A,,(p) and A(pl + a) are, respectively, elements of the adjunct matrix and the determinant of the matrix 
pl +a. We elements Bj, are found from the formula: 


al (T'—a) 
“ji, = x 


We consider the general case, when the characteristic equation A(pI + a) = 0 contains multiple roots, In 
this case, 


N N 
A(pl +a) P+)", 


vel vel 


The operator functions in (24) can be represented in the form of sums of elementary fractions, 1.e., 


(25) 
Bolsa (Pp +2)" 


where 


The adjunct matrix A(p) and its elements possess the following properties (Cf. [12], [13], pp. 74-80): 


Aj (— dv) = ASR (— 4) =... = AN (— = 0, 


(26) 
Atk (P) = An (Pp). 
On the basis of (26), we get 
dy 
Kk — Kiko... = Kik =0. 
By substituting (27) in (25), we get 


In Equation (28), we now make the transition from the transform to the original time domain, taking (23) 
and (24) into consideration, As a result, we obtain 


N 


vel 


Consequently, 


N 
— Gj, k= 1, 2 (29) 


vel 


In order to put the matrix e~*T® in the form of (29), we used the well-known operator transformation which 
underlies the transformation to canonical variables and normal coordinates [11]. 


By substituting the value of e *?" in Equation (22) and presenting it in expanded form, we get 


al 
| 
x 
in 
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N r 
KU m 10) + 


vel mel 


2K yk, j1,—A yT (n—k—1) 


vet k=o 


In the case when the step function Aafn] has a constant value, i.e,, A [n] = Ad, the expression for Uj{n] 
takes the form: 


N r 
= KU m 10] + 


vel m=1 
i—e 
+ a Ky pe ey 


Expressions for Uj[n] are easily found in other cases, if Ao{n] is given by elementary functions, 


7. Discrete Equations of the Auto-Range-Finder 


For an auto-range-finder considered as a closed system of automatic control, the discrete equations for 
the general case are written in the form: 


equation of the time discriminator and operational amplifier 


r n—1 
= Brn Um {0] (nT — AT —7) (30) 
k=0 


equation of the regulated delay circuit 


A$ [n] = [n]; (31) 
equation for the auto-range-finder tracking error 
Aa [n] + Ad [n] = Ad, [n]. (32) 


To obtain Equation (30), we employed Expressions (22), (20) and (23). Equations (31) and (32) were set 
up in correspondence with (4) and (5). 


8. Operational Representation of the Auto-Range-Finder Discrete Equations 


To solve the system of discrete equations obtained in paragraph 7, we use the method, presented in the book 
of Ia. Z. Tsypkin [6], of the operational calculus of step functions, We shall denote the discrete Laplace transform 
of a step function f[n] in the following way: 


L* = p=e+io, e>0, i= 
n=0 


We now apply the discrete Laplace transform to Equation (21). Using the prediction theorem (Cf. [6] pp. 
21-23), we obtain 


(33) 
L* (U [n}) = 10] + (Ax 


Corresponding to (16) and (20), if Inequalities (17) hold, the matrix e~*T can be given by the approximate 
expression e~*T w1—aT. By substituting this expression in (33) we get 


2K 
[PI +a] L* U = [0] + (Aan). (34) 
eP — 1 
Here, P=—7—. 
We find from (34) that 
(P) {0} 
Ant 


Corresponding to (31) and (32), we also have: 


L* {A$ [n}} = oL* {U,[n}}, 
{Aa + L* {AS [n]} = L* {Ad [n)}. (36) 


By solving Equations (35) and (36) with respect to L* {Aaf{n]}, we get 


[aa + (P)] (Aa = 


= Aa(P) L* [n]} — 5 Arm (P) Um [0], (a7) 
m=1 


Aa(P) = A(PI +a). 


We now consider a continuous system with constant parameters, the processes in which are described by the 
equations: 


[D + (t) = Aa (), 


> D + Ux (t) = 0, (38) 
k=1 


AS(t) =U, (t); Aa (t) + Ad (t) AS, (2). 


Applying the Laplace transform to Equations (38), and solving them for L[Ac(t)], we obtain 


Aa (p) =A (pl + a). 


(39) 


By comparing Equations (37) and (39), we find a relationship subsisting between them, thanks to which we 
may simplify the task of finding the representation of the step function Aaf[n]. 


9. Transition to the Discrete Equations from the Continuous Ones 


We set 


(40) 
L*{A9, [n]) = Aa(P)@,(P) => 
In addition, we introduce the notation 


20K. 
X(P) = Aa(P) + Ay (P). (41) 
The functions X(P) and Y(P) are polynomials with constant coefficients: 


X (P) = + 


We shall assume that the characteristic equation X(P) = 0 contains only simple roots, In this case, the 
function X(P) takes the form 


X(P) = 


Ver] 


where the yy are all different values, 
Solving Equations (35) and (36) for L* {A a{n]}} , and taking (40) and (41) into account, we get 


Anm(P) 


We present the fraction in this equality in the form of a sum of elementary fractions: 


The constant coefficients Gy, Gy and Sy are found by well-known formulae (Cf., for example, [13] chapter 
V1). 


cone (43) into account, we made the transition in Equation (42) from the transformed representation to 
the original,” As the result, we get 


Ma [n] = + G,+ 5.(1 + pT)" 30.101 (4 + 


v=1 vel 


X' (P) = 4 X(P). 


Trables of inverse transforms, given on pages 211-217 of la. Z. Tsypkin's book [6], can be used in making the 
transition from the transformed equations to the original ones. 


: 
| 
nm 
(43) 
Here, 
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In the general case, the roots py can be either real or complex. If the following inequalities hold for 


then? 
[RepT |< 1, |Imp,T|<4, (44) 
then 1+ py T 
Consequently, 
r 
Aa [n] = + G+ Sve" NUS (45) 
| 


In the case considered, the step function Aca[n] can be obtained from the continuous function Ao) in the 
following way: 


Aa [n] = [Ag ()lenr (2 =0, 4, 2,...,00)) 


Aa (t) = Got + G, + —o 
vel vel 


The function Aogt), in its turn, is found from the solution of a system of equations with constant coefficients 
of the form of (38), if L[®(t)] = OQ(p). In other words, if Conditions (44) hold, the investigation of processes in the 
auto-range-finder can be replaced by an investigation of the continuous processes in a system with constant para- 
meters. This turns out to be valid, not only with respect to the variable Aafn], but also for the other variables, 
Uy[n] and A®[n]. One can convince oneself of this by carrying out similar computations for the other variables 
also, computations which were omitted here due to their cumbersomeness, 


Similar results can also be obtained for the case of multiple roots, if one uses the form of solving systems 
of linear differential equations obtained in the work of B, V. Bulgakov [11]. 
SUMMARY 


The auto-range-finder equations, containing variable parameters varying step-wise, transform, in the general 
case, to a system of discrete equations with constant coefficients (30), (31) and (32), the solution and investigation 
of which pose no particular difficulty. Conditions were found for which the analysis of the auto-range-finder func- 
tioning reduces to the analysis of processes in some equivalent continuous system. 
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ON SERVO SYSTEMS CONTAINING TWO PULSE ELEMENTS WITH 
UNEQUAL REPETITION PERIODS 


Fan Chun-Wui 
(Moscow) 


Pulse servo systems containing two pulse elements with unequal repetition 
periods are considered. Their equations and transfer functions are found, An investi- 
gation is carried out on the influence of the inequality in repetition periods of the two 
pulse elements on system stability. It is shown by an example that, in the general 
case, Ty = Ty (equal repetition periods of the two pulse elements) is not the most ef- 
ficient regimen for the system from the point of view of increased stability, 


Previously, systems were investigated which contained one or several pulse elements with identical repe- 
tition periods [3-5, 7-9] or with repetition periods which were simple multiples of one another, Systems containing 
several pulse elements with arbitrary repetition periods remained uninvestigated, despite the fact that this problem 
has essential theoretical and practical value. 


1. Transfer Function of a System Containing Two Series-Connected Pulse Elements 
With Unequal Repetition Periods 


Figure 1 shows an open-loop system with one pulse element, We denote by H(S) the transfer function of the 
linear portion, and by h(t) = L™![H(S)] the impulsive response function of the system, If it is assumed that the 
porosity of the pulses yp << 1, then the transfer function of the pulse system, with the condition that h(0) = 0, will 
have the form [6]: 


H* (STg, = H(S + 


l=—o 


where Ty is the repetition period of the pulse element, w= 2a/T, and the parameter € varies within the interval 
O0se<1, 


It is clear from (1) that the transfer function of the system is 


H for T = 
7] 

> H(S+ jlo.) for e=OnT+0, 

© 

t=—o 
h(o 


** That is, without a pulse element, 


where H*(STp, 0) is the discrete Laplace transform of the function h(t) = L*[H(S)j, {.e., 


(ST,, 0) = D{h(t)} = h(mT,) e-St, 


m=0 


We now -onsider the case of two pulse elements, connected in series, with corresponding linear portions 
(Fig. 2). The sollowing relationships hold: 


ae X (S) = C€(S) H*(ST;, ©), C(S) = E*(ST,, 0)G(S), 


(ST, 0) = H(S + 


and Ty and Tg are the repetition periods of the first and second pulse elements. 


| FE ms) PE Gis) PE Ms) 
Fig. 1. Fig. 2. 


From the previous relationships, it follows that 
X (S) = E* (ST,)G(S) (ST, = 


= E* G(S) A(S + jleg) (2) 
l=—oo 


If the discrete Laplace transform is applied to the function x(t) = L™*[X(S)] at the points mT}, it is neces- 
sary to replace €T, by the values mT, in Equation (2). 


Since E*(ST,) is a periodic function, with period w, = 2n/T,, 


X*(ST;) = X(S+ = 


= jme,) H (S + + jmo,) 


The transfer function of the system in the open-looped state is given by the expression 


G(S + jme,) H (S + jlo, + jmey) 


(3) 


For the system in the closed-loop state (Fig. 3), we have the equality: 
E* = F* (ST,) — X* (ST)). (4) 


The transfer functions of the closed-loop system, determined from (3) and (4), have the form: 


(S, T,, 7.) = X°(ST:)_ (S, T1, Ts) 


(5) 


F’(ST;) i + KW (Ss, Ti, 


We consider the following particular cases, 
1, For Ty = Tg = To, Expression (2) takes the form: 


RW" (ST,) G(S+ jmoy) H(S + (ST_) (ST) 


2. For Ty = Ty and Ty = 0 (i.e., without a second pulse element), it follows from (2) that 


KW” (ST.) = 7 G (S + H (S + 


It is clear, from the argument presented above, that 
Formula (3) expresses the general transfer function of an open- 


19) BE 
ea eae looped pulse system, Since the linear portions, with transfer 


functions G(S) and H(S), are generally low-frequency filters then, 
for sufficiently large values of wy and ws, it is possible to write 


Fig. 3, 


> 2 G(S + (S + + 


M=—co l=—co 
—j on 
G(S)H (8) +G(S+ jo) |H(S + joy) + HIS + 
ja 
"| + G(s — jo) — jo) + 
on 


(6) 


This formula permits the frequency characteristic to be constructed approximately, after which the stability 
and quality of the pulse system can be investigated for arbitrary values of Ty and Tx. 


2. Equation of a Discrete Correcting Device Consisting of Two Pulse Elements. 
The Repetition Period of the Second is a Divisor of the Repetition Period of the 
First 


We assume that the repetition period of the second pulse element Tg is an integral divisor of the first, 
i.e., Ty = Ty and T, = Ty/n, where n is an integer, In this case, (3) takes the following form: 


RW*(S, + H[S + jos(t-+ (1) 
m=—co 
In this case, the linear portions of the pulse system, with transfer functions G(S) and H(S), can be inter~ 
changed without changing the over-all system transfer function, To convince ourselves of this, we prove the 


validity of the following identity: 


(8) 


+ imo G[S + jor(t + 


To prove this last identity, we assume that (2 + m/n)w, = Muy and mu =(L + M/n)wg. It follows from this 
that | =—L and m = nL + M. 


By substituting the values of 2 and m just obtained in (7), we find that 


nL+M=-—co L=—co 


+ i(t+ coy] H (S + 
From other considerations, we obtain (8). 
It follows from Fig. 2 that 


C(S) = E*(ST,)G(S), X (S) = C*(ST,) H(S). 
It is known that 


C(S + jlo) = 


+ jl) G(S + jlo) = 
l——oo 


=F + + jlo) G(S + jlo) = 


m=]—oo 


LS + + jlon) = 
m=—a 
+ ima) G(S + = 


= E*(ST,)G* (ST,). 


By substituting the expression for G* (ST ) in (9), we find: 


X (8) = E*(ST,)G" H(S),- 
X*(S7,) (S + jma) = 


1 


=F (87) 2% + + + 


Expression (7) can be rewritten in the form 


X*(ST,) = (ST) G(S + jmoy) H[S + jo (t+ 


By combining the last two expressions for X*(ST 4), we obtain Identity (8), q.e.d. 


We now assume that H(t) = L“{G(S)H* (ST,)/n)] and (t) = L™{H(S)G*(ST)/n)]}. Identity (8) shows that the 
discrete Laplace transforms of these two functions are equal, i.e., 


KW* (ST,) = D (t)} = Din}, 


or 
T 
(ST) = D{L~ 1G(S) = D{L 
The last two expressions for K*(ST,) corresponds to 
| Fig. 4,a and Fig. 4,b. If, in these diagrams, we introduce 
= PE | GIs) a third pulse element with period T,/n before the element 
f with period Ty (Fig. 4,c), the systems considered do not 
| change their properties and, consequently, the last expres- 
| sion can be rewritten in the following form: 
[PE | 
| KW* (ST) = D [G* (ST, /n) H* 
| where we denote by D,, “! the operator for the inverse dis- 
ea) (PE Rute PE PE wis crete Laplace transform, and (ST,/n) and H*(ST,/n) 
} are the discrete Laplace transforms of the functions g(t) = 


= L~*[G(S)] and h(t) = L™! (H(S)] with period T,/n, i.e., 


ST 


G(s 72) = Dee 


72) = Dy {h(t)} = 


If the transfer functions of the linear portions are known, 


then KW*(ST,) can be computed from the formula (Cf. Appendix): 


bj A, — 


jaye n A,T, —b;, A,T, 


where T is a constant value of lag, and 44, 4, and mg are as shown on Fig. 8, 
In the case when some poles of G(S) and H(S) coincide, i.e., for 


901 


| 

| 

| 

Fig. 4, 

| 

| 

d 

= 


d d 
H (S) = 


Kw*(ST,) can be given as the sum of two terms KWw(ST,) = KWy*(ST,) + KW,°(ST,) the second of which, KW,°(ST,), 
(the sum of products of terms with nonidentical poles), is determined from Formula (10) and the first of which 
Kw¢ (ST) (the sum of products of terms with identical poles), is determined from the following formula (Cf. 


Appendix): 


KW; (ST) = 
nr e(l—m,)8T, [ + As — )| (11) 


In this case, the transfer function of the closed-loop system has the form: 


X*(ST,) [o (s H* (s Kw*(sT,) 


~ T+ RW’ (ST,) (73) 


It is possible, using (10)-(12), to carry out an analysis of the quality of the choice of correcting devices 
for such systems, If we are interested in the output quantity, not only at the instant t = mT», but also at the 
instant t = mT,/n (mand n integers), then (12) can be written in the following form: 


(s 72 F*(ST,) = 
Kw* (s =) F*(ST,) 


This expression coincides with the result given in work [1]. 
Using the notation X*(ST,/n)/F*(ST,) = ;°(ST)/n), we write (13) in the form 


Rw’ 


7) = 


where KW®(ST,/n) is the desired function. 
The solution of this equation (as shown in[1]) has the form: 


Whence, 
o° (s To 


The relation (14) determines the transfer function correction for the transfer function object H (S) and 
functions sis. and (ST,). 


3, Investigation of the Influence on System Stability of the Inequality of the 
Repetition Periods of the Pulse Elements 


As an example of the use of a discrete correcting device of two pulse elements (the repetition period of the 
second being an integral divisor of that of the first), we consider the system whose block schematic is given in 
Fig. 5. We assume that the system possesses the following 
properties: 1) it has first-order astatism and 2) it provides 
Als) pas PE the minimum possible regulation time, It is required to 
determine the form of the correcting device, 


We consider the following cases, 


Fig. 5. a) n=2,T, =1 second, We determine H*(ST,/2) 
in the ordinary way: 


e* —0,37 


It is known [1] that for a system with astatism of the y ‘th order, the transfer functions @*(ST,) and ,°(ST»/n) 
have the following forms [1]; 


—8T.\w 4° 
1—0*(ST,) = “=! YA" (STs) 


where A*(ST,) and B*(ST,/n) are polynomials in e°T and eSTo/ " and the cy are real coefficients, In the given 
case, n = 2 and y = 1. Condition 2) requires that Ay(ST,) = 1 and cy, = 0, In this case, we obtain 


(16) 
(ST,) = 


Since = 0" (ST), ten BY ($28) 


= 
q 
q 
| 
*) 
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By substituting the corresponding terms, we get the schedule for computing the correcting device: 


8T, 17 
(5 (17) 


0; (842) = 


(18) 


b) n=1, Ty = 1 second. In this case, (14) takes the following form: 


1 *(ST,) 
G* = (ST, i— 


The functions *(ST,) and H®(ST)) are easily determined: 


0.86 


=e", H* (ST) = 


From the relationships given above, we obtain: 


1.2e5T+ — 0.16 
G* = (19) 
c) n=1, Ty = 0,5 seconds, Analogously to the previous case, we have 
1.6e —0.58 
G (s = (20) 
e*—i1 


We remark that (20) coincides with (18). 


Curve 1 of Fig. 6 is the output function of the system X(mT/2), 
where mis an integer, for unit external excitation, for the cases a) and 
a c), and curve 2 is for case b). 


It is clear, from the example given, that cases a) and c), trom the 


bull will be simpler for case a) than for case c). 


— a static system of automatic control. We assume that, in this system, 
K i—e * 
G(S) = and = 


We now find the maximum attainable value of system gain, ee for Ty = Ty and T; = T,/n (n an integer). 


By means of Formula (10), the transfer function of the open-looped system can be determined: 


K (c,e8"* + 03) 


KW* (ST) = 


point of view of rapidity of action, are equivalent (Cf. curve 1 on Fig. 6). 
Initial alignment will have the same transfer function for cases a) and c), 


12@9%465 678 9 0m In the second example, we investigate the influence of the inequality 
of the pulse elements" repetition periods on system stability. We consider 


| 

ST, 

—1 
MT) 

904 


where 


Ts 
pls 


To simplify the computations, we assume that T = To, i.e., mg = 1, Ay = 0 and 4, =1/n. By substituting 
these values in the expression for cg, we get 


Ts 


Then, the transfer function of the open-looped system takes the form: 


Key 


and the characteristic equation of the system is 


By employing the analog of the Routh-Hurwitz criterion, we determine the following stability conditions; 


1. (4—e*™) (1 + > 0, 
3 (4+ (1 >0. 


Since Cc, is always positive for any positive values of a and b, Conditions 1 and 3 always hold, Therefore, 


; the condition for system stability will be 
—(a+b)T, 
Rmax< 
ality 
ler 


_With given values of a, b and Ty it is possible, from this last inequality, to construct the functional relation- 
ship, Kmax = f(n) (where n is an integer). For intermediate points, it is necessary to use Formula (6), since 
inequality (23) will not be valid at these points. The following table gives values for c, and Kmax- computed from 
Formulae (21) and (23), for various values of n, for a = 4, b = 2 and Ty = 1 second. 


n i 2 3 4 5 6 7 8 9 10 oo 
Cy 0.87 | 0.32 | 0.22 | 0.19 | 0.18 | 0.16 | 0.155 0.150) 0.145 | 0.140) 0.117 
1.16 | 3.18 | 4.50 | 5.16 | 5.65 | 6.16 |6.44 [6.76 | 6.90 | 7.10 | 8.55 


The curve for K,,,,, = f(n) is given in Fig. 7. 


le le 4 
~a b 
q 
q 
| 
d 
j 
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It is clear from these examples that, with increasing n, f.e., 
with a decreasing repetition period of the second pulse element, the 
stability margin of the system increases in comparison with the case 
Ty 


ethod f onstructi 
Response in Pulse Systems Containing Two Pulse 


Elements 


The transient response can be determined analytically by 
means of the formulae for inverting the Laplace transform, or con- 
$6676 9 tn structed by graphic or tabular fitting. 


Fig. 7. We investigate both these methods, 
1, We use the formula for inverting the discrete Laplace transform [2]: 


—— 


r, 


Here, n is an integer and the integration is performed along the contour Ty, bounded by the lines corres- 
ponding to c and jx/u. The function X*(ST, « ) occurring in (25) can be found by applying the discrete Laplace 
transform according to Formula (2): 


X(ST = PR H(S + + jmoy) x 


m=—o 


wos 


From these last two decompositions we obtain: 


x 
on jim ox 
TT, + H(S + (36) 


x em (n+*)T, d 


In the general case it is possible, in principle, to obtain the transient response by using the last decompo- 
sition, However, the integration can be simplified if approximate Formula (6) is used. 


2. It is clear from Fig. 2 that 


X(S) =C*(ST,) H(S) = HS, 
l=—co 
(S) = E*(ST:)G(S)= (5), 


m=—coO 


(/T 2) = c(t) = [Cc S)] 


e(mTs) = €(t) = SM 
By using the formula for the inverse Laplace transformation, we obtain 


z(t) ¢ (ITs) h(t 


c(t) = e(mT,) g(t — mT,). 


From these last two expressions we get 
nT;/T, 
z(t)= > g — h(t — IT). 


l=0 m=0 
The last expression can be given in the following expanded form: 


z(t) = h(t) [e (0) g(0)] + 
+ h(t —T,) [e (0) g (72) + (71) g + e (273) g (72 — 273) «J+ 
+ h(t— [e(0) g (272) + €(7 3) g(2T.— 7s) + g(2T,—2T7)) +...J+ (21) 


+ h(t—KT,) Je(0)g (KT) + + 


Such is the method of attacking ordinary pulse systems expounded in [2, 4]. 
If the external forcing function f(t) is given, and g(t) and h(t) are known then, by using Formula (27), it 
is possible to construct the system's transient response step by step. 


Appendix 
We assume that G(S) and H(S) have the form: 


d 
G(s) = H(S)= 
where T is a constant value of the object lag. Then, we can find KW*(ST,) = D {p7tc* (ST»/njH* (ST»/n)]} 


in the general form, Figure 8 shows a typical transientresponse h(t). We assume that T = (mg —44) (To/n), 
where my is an integer and 0 < Ay < 1, We find G*(ST»/n) and H*(ST,/n): 


e 

4 

{ 

e *—e 

Ts 

8S — 

e —e 
{ 


We use the notation: 


n 
—s 2*-m, a ats sis 
VA 
/ 
Fig, 8. 
Knowing KW®(ST»/n), we can find the corresponding function: 


Since D {D;" [KW*(ST,/n)]} = KW*(ST,), we can from this, obtain the generalized formula for KW*(ST,) 
(Cf. Fig. 4): 


—b,A, Te 
n ,—m,ST, 


where m, =(m,—1)/n + Ag, mg is an integer and 0 <A, < 1 (Fig. 8). 
In the case T = 0 (i.e., my = mg = 0, Ay = 1/n), the last decomposition assumes the following form: 


e —e e —e 


We now consider the following particular cases of the last expression. 
1, For n = @ (i.e., Tz = 0), Decomposition (29) takes the following form: 


els 


») 


2. For n = 1 (i,e,, Ty = Tg Tg), it will be 


Expressions (29) and (28) are also valid in the case when there is one pole at the origin of coordinates in 
G(S) or H(S). If some poles of G(S) and H(S) coincide, i.e., 


G(S) = 


+ 


then in this case KW*(ST)/n) can be given, as previously, as a sum of two terms; 


agit a, 282! — bj 
[a> “e e od,e "e e aly 
Ts T Ts 
+ 


where KW,°(ST,/n) is a sum of products of terms with nonidentical poles, The form of each of these products 
coincides with Expression (28), so that is necessary to consider only KWy (ST,/n): 


Ts Te 
os A, 2? 


We find from this last expression that 


1(S)= Ya ‘ 


Knowing the expression for KWq(S), we can easily find 


7 J 
4: — T, T. 


d 
(ST») = Dy 


—a, 
e +) 


| 

and 

a 

7 
| 

| 

i 
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For T = 0, the last decomposition takes the following form: 


4 ST. —a,T. 
e —eé 


It is easily seen that this last expression also holds in the case when one of the roots a) is zero, 


SUMMARY 


1, The transfer function of a pulse system containing two pulse elements with unequal repetition periods 
was found, 


2. By using the given simple conditions for a correct choice of the ratio of Ty and Ty, one can increase 
the stability margin of the system, In general, the condition Ty = Ty is not the best regimen for system functioning 
from the point of view of increased stability margin. 


The author wishes to express his gratitude to Professor la, Z. Tsypkin for pointing out the direction that the 
investigation of this subject should take and for his valuable advice, 
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DETERMINATION OF PERIODIC BEHAVIOR OF AUTOMATIC CONTROL 
SYSTEMS CONTAINING A NONLINEAR ELEMENT WITH 
BROKEN-LINE CHARACTERISTIC 


L. A. Gusev 


(Moscow) 


A method is considered of determining the periodic behavior in such 
automatic control systems which differ from the linear ones by the presence 
of an arbitrary characteristic of broken-line type. The periodic behavior is 
fully described by complete Fourier series, without neglecting higher 
harmonics, The problem is reduced to the solving of h simultaneous trans- 
cendental equations (equations of periods) which determine times of motion 
in each part of the nonlinear characteristic within the period limits, Also 
some problems connected with the use of computers for solving period equa- 
tions are considered at the end of the paper. 


Self-oscillatory behavior is characteristic of some automatic control systems, Consequently it is often 
necessary to determine the period and the form of self-oscillations from the known equations of all elements of 
the control system. In systems having vibrators, or in systems subjected to periodical perturbations, there arises 
the problem of determining the period and the form of the forced oscillations innonlinear systems, An approxi- 
mate answer to these questions, for example, by harmonic balance method, often proves unsatisfactory as the 
automatic control systems do not by any means satisfy everywhere the conditions under which such methods are 
applicable. In such cases the periodic behavior must be calculated exactly, and this always involves considerable 
difficulties, Even when one succeeds in separating one “principal” nonlinearity in the system, and in reducing 
the problem to the study of an automatic control system which differs from a linear one only by the presence of 
one nonlinear element, even then the exact solution can only be obtained on the assumption that the character- 
istic of the equation is of a broken-line type. The problem is reduced to the solving of period equations, that 
is, of final — even if transcendental — equations whose roots, while satisfying some additional conditions, deter- 
mine directly the required periodic behavior. There are two methods described in the literature which permit 
exact determination of the periodic behavior in this broken-line case, that is, permit the construction of period 
equations, 


1, “Make to fit” method, In applying this method it is necessary to integrate out all equations of the linear 
systems into which the nonlinear system under consideration has been decomposed, and the periodic solutions are 
constructed by using these integrals. 


2, Finding periodic solutions in complete Fourier-series form (not neglecting higher harmonics), In this 
case it is not necessary to integrate the corresponding linear equations, and the period equations are directly 
expressed by coefficients of original simultaneous equations which describe the motion of the automatic control 
system under consideration, The problem which has been solved by this method was a particular one when the 
characteristic of the nonlinear element consisted of finite straight lines parallel to two given straight lines only 
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The method is applied to more general systems in this paper, that is,to systems with an arbitrary broken-line 
characteristic (the number of different directions of portions of the characteristics is now arbitrary). As distinct 
from [1], the affine transformation of the given characteristic into a characteristic consisting of straight lines 
parallel to coordinate axes is not applied here. This brings about a considerable simplification in the equations 
which serve to determine the times of motion of a representative point in each separate portion of the character- 
istic, 


If computors are available it is convenient to apply the proposed algorithm when there are a few component 
portions of the characteristic and the linear part of the system is of a high order. In a system with a low-order 
linear part it is obviously simpler to have a model of the problem and “to get the feel" of the phase space, and 
thus detect the stable limit cycles. 


The stability of periodic behavior is not discussed in this paper. These problems are considered in [2]. 


1. Introductory Relations 


Our problem can be mathematically formulated as follows, We have simultaneous differential equations 


dz 
= > xy, Y= (i = 1,2, ...,n), (1) 
j=1 


where Oy and K, are real, and y = f(x) is any broken-line function, whether continuous or discontinuous, 


The transition from one portion of the characteristic y = f(xy) to another, is said to be taking place normally 
in the sense of [3] when sliding or ambiguous transitions are absent, The type of periodical solution, that is, the 
sequence of successive straight lines which form the characteristic y = f(x), is given for the duration of periodic 
motion, The problem is then reduced to the solving of the final period equations which determine the times 
taken by each of these portions. 


Eliminating from simultaneous Equations (1) all x; except X;, we obtain equation 


D (p) = K (p) f (2) (2) 


which we shall call the derived equation (to simplify notation, subscript 1 of Xq shall be omitted from now on). 
In addition, for every moment of time t_,during which a break occurs in any function X(t) or function y = f(x) or 
in any derivatives up to order (n —1) inclusive, we obtain gap relations (see [4]): 


= 
+On—1 bo = +--+ + No- 
In (2) and (3) 
D (p) = + + an 
K (p) = bop" + + On 
are linear differential operators,* and 
= 2" (tg + 0) —2 (tg —0), 
Tr = y™ (tg + 0) —y™ (tg — 0) (4) 


* The degree of K(p) is less than that of D(p) but we regard them as equal by assuming that the highest coeffi- 
cients in bg, by,.. ., bg are equal to zero, 
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are the heights of gaps of the k-th derivative (k = 0, 1,... , m1) of functions X(t) and y(t), respectively, at 


the time t = tgs 


Further, quantities €, and ny, are considered bounded, that is, no gaps of the second type occur either in 
the functions or in their derivatives, 


When the equation of the investigated straight-line piece of the characteristic with subscript q is written 
in the form y = kqx + Ag and then this y is substituted back into (2), we obtain 


Lq(p) r= Cy (q = 1, 2, A), 


where the notation Lq(p) = D(p)—kqK(p) = gop" + + By = byAgs BE = 


The pattern of the required periodic behavior determines the sequence of the successive differential equations 
of the type of (5). Some equations in this sequence, though corresponding to different subscripts, may be identical, 
due to the fact that during the periodic motion one and the same portions of the characteristic can occur more than 
once, 


(5) 


Therefore to solve our problem it is necessary to find a periodic function x(t) with period T = &, and satisfying, 
successively, the equations 


Ly (p) z(t) =e when hy, 
Ly (p) x(t) = cg whent, <t< te, (6) 


Ln (p) (t) = cy whenty_y<t <tr. 


Then the gaps of the functions x(t), y =f[X(t)] and of their n —1 derivatives should satisfy the gap relations 
(3) at each point tg(q = 1, 2,...,h). The numbers ty, tg... , t, are unknown and must also be determined, 


2. Transforming Equations andGap Relations 
We shall find the periodic solution of Eq. (6) and (3) in the form 


h 
x(t) = x(t) + pull (1) 
q=1 


h 


where z(t) = bs Zq(t), 1(t) is the Heaviside unit function, pg = 4 » and x,(t) is an auxiliary periodic 
q=1 n 


function, with period T, and satisfying the following conditions; 


%q(t)=0 when OCt*< ty, 
Lq(p)%q(t)= 0 when ly 
aq(t)=0 when ty <t< the (8) 


The sum on the right-hand side of (7) represents a step function of t, taking the values Pq in the intervals 
t <t<t 
q~1 


For a particular interval with subscript q the transition from x(t) to x(t) in accordance with (7) is in the form 
Z(t)=2(t)t pg (9) 
Therefore the q-th equation of the System (6), valid in the same interval, takes the form 


Lq(p) (t) = Lq(p) (x(t) + pq] = Le (t) + = Cg. 
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But pq = cq/gh and when t,_y < t< tq, x(t) is determined by the homogeneous equation 


z(t) =90 (qg=1,2,...,h). (10) 
Introducing the notation 


(ty +0) 28, 0) = 
(ta —0) = (tg —0) = (11) 


and taking into consideration the equation of the characteristic y = f (%) and Relations (7), we are able to trans- 
form the gap relations, As a result we obtain a system of n linear equations (deduced in Appendix 1): 
= Ug 2, + — 
q2 (12) 


Here we have introduced the notation = All the other and vg are 
given by recurrence formulas 


+1 ok. Pa 
qk 


where Pox = (bx — — pq (de — 87%). 


We now remember that by and gi! are the coefficients of the polynomials K(p) and Lg(p). As all the 
functions considered here and below are periodic, we are allowed to change subscript h + 1 simply into 1, 


We shall call System (12) the cross-over condition, from the q-th to the (q + 1)th piece of the character- 


istic, 
Introducing column matrix 
x(t) 
|, 
(14) 
X (tg 0) => Xa = 
we write (12) in matrix form: 
Xg=UXa—Ve =1,2,...,h). (15) 


In (15) Ug is a triangular truncated-rows matrix: . 


are able to evaluate det U, by making use of the footnote on p. 912 From either Ugg = +8709 ay/ay = 1, 
or by = 0 we have det Ug = 1. Therefore the matrix Ug is nondegenerate, k, = @ could be an exceptional case, 
but then as long as the motion is in the portion corresponding to k_, X = const, that is, it is known, 
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‘Ugo 
Uqi Ugo 
U, => Ug 


oo 


and is the column matrix 


Va 
Ve = . 


Vq, n—1 


Thus, the previously formulated question has been reduced to the determination of a discontinuous periodic 
function x(t) satisfying, successively, Equations (10) and cross-over Conditions (12) or (15) at every point of dis- 
continuity t = tg. 


3. Construction of Auxiliary Function Xq(t) 


We are introducing periodic functions Ryq(t) (k = 0; 1,...,m-~1), period T = t,, satisfying the following 
conditions: 


1. Function Ryg(t) is a solution of the differential equation Lo(P)Ryg(t) = 0 for all t except t= ¥T (ves 


2. All derivatives Re) when 0 <j < kare continuous, 
3. k-th order derivative R(t) is discontinuous at points ¥T; the magnitude of its gap is unity, 
4, All derivatives Ry when k< j = n~—1 are again continuous,* 


The construction of functions Rxq(t) is shown in Appendix II, and it is also shown there that their Fourier- 
series expansion is of the form 


(16) 


Ig(s) = gt_, + + 


Making use of the obtained functions Ryg(t) we construct the periodic function x,(t) defined by the expres- 


sion 
n—1 n—1 
tq(t)= >) Rig (t — tyr) — Rag (t — te) Br, (17) 
k=0 k=0 
where a) and 6, are arbitrary real numbers and 0 _, < T. 


l, —— 
* The discontinuous function Rn is not differentiable at the points of discontinuity, and the continuity of its 
derivatives is understood here in the sense of formal continuity, that is, the left-hand and right-hand side deriva~ 
tives are equal to one another at every point of discontinuity vT, 


Rig (t) = 
where s==irw (w= Fr), and 
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On the basis of the properties of functions Ryg(t) one is able to draw the following conclusions regarding 


the properties of Xq(t): 
4) Xq(t) is a solution of the differential equation Lq(p) Xq(t) = 0 for all t with the exception of points t 4+ 


+ VT and tq + UT (integral ¥) and, therefore, the function and its first n—1 derivatives are continuous everywhere 
except at these points; 


b) within period T the gap of the k-th derivative of the function Xq(t) at the point tq—4 18 equal to ay: 
(tg—1 + O) — (tq—1 — O) = a. 
Similarly, at the point tg we have 
(tq + O)— (tq — O) = — Br. 


Having evaluated X(t), Xq(t),..., xn —) (1) from (17) (the evaluation method of the derivatives of function 
Ryq(t) is described below), and having formed the column matrices 


t 
(n—1) (t) 


we are able to write the obtained equalities in matrix form: 
Xq(t) = Mg (t — tg) A— My (t — ty) B, (18) 
where Mq(t) is the matrix 


Ro (t) Rig (t) a(t) 


Rit) 


Now it is required that the function X q(t) defined by (17) be identically zero in the intervals (0, ty _4) and 
(tqr T) that is, that it be the trivial solution of equation Lq(P) Xq(t) = 0 in these intervals, It is obviously sufficient 
tor it to satisfy conditions x,(t) = Xq(t) 2...EKy (t= 0, or in other words, xf = 0 for any t from these inter- 
vals, In particular, having taken t in the interval (tg, T), we obtain x_(t) = 0 not only inside interval (tg, T) but 
also inside interval (0, tg _s) because function Xq(t), together with its n—1 derivatives, is continuous at points 0 


and T, according to property (a). 


We select a point tq+ 0 in the interval (tg T). At this point the requirement Xq(tg + 0) = 0 can be written 
in the form 


Mg (tq — A— Mq(4+0)B=0. (19) 
in agreement with (18). 
The values of A and B satisfying (19) will now be considered. We have then 
Xq = Xq (tg 0) — 0. (20) 
From property (b) of function Xq(t) we have 
Xq(tq—1 + 0) — Xo =A, (tq —9) = —B. (21) 
It follows from (20) and (21) 
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Xq(tg1+0)=A, Xq(t*—O) =B. 

The choice of A and B satisfying (19) ensures that they are equal to column matrices of initial and final 
values of the function Xq(t), (an integral curve of the equation), and of its n ~1 derivatives at both ends of the 
interval (tq 4+ tq). Therefore one can substitute for column matrices A and B the column matrices Xqltg ~4* 
and Xq(tq —0) in (19) and (18), Having performed the substitution, we obtain 


M gq (tq — Xq + O) Mg (+0) Xq (tg —O) = 0, 


(22) 

Xo lt) — + 0)— Mg (t — tq) Xq(tg—9), (23) 
&q(t) = Rug (t — tes) + 0) — Rag (t — tq) 20 (tg — 0). 

k=0 k=0 (24) 


Formulas (23) and (24) together with (22) determine the periodic function xq(t), period T = q,, satisfying 
Conditions (8). 


4. Period Equations 


We are similarly constructing functions Xq(t) in each interval (tq t) (a= 1, 2,...+,h), The sum x(t)= 
[and correspondingly X(t) = Xq(t)] complies with (10). 


Now it is necessary that the gap Conditions (15) be fulfilled at each point tq. Taking first into considera- 
tion the obvious equalities [see notations (11) and (14)] 


X (tg —0) = Xq = Xq(tg—9O), 
X (tq + 0) = = (tg +O), 
X + 0) = = Xq (tgs + O) 


and then making use of (15) we can put (22) in the following form: 
Mg (tq — ta—1) X q—1.2 — (+ O) gX qn — Ma (+0) = 9. (25) 


Relation (25), when the value of x(t) and its derivatives at the point toast 0 is known, permits to determine 
the values of x(t) together with its derivatives at the point tat 0. 


Taking into account the gap conditions, formula (23) becomes 
Xq(t) = Mg (t —tq-1) Xq—1.2 — — te) UgX en — Mg (t — tq) (26) 


By adding to the system of Equations (25) the condition Xp. = Xgg, which should be true in view of the 
periodicity of the required solution x(t), a closed system of equations is obtained: 


M, Xoa — Ni (+ O) = S, (+ 0), 
(te — X12 — No (+ O) = 0), 


Mp (tn — tna) Xn—-1, 2 — Na(+ O)Xna = Sp (+0), 
Xne = Xo, 
where 
M,(t)Ug= N,(t), M,(t)Vq = S(t). (38) 


If ty, ty, » t were known, one could find all = Xqa(ty, tg, - + » %,)from System (27) and construct 
all functions , that is find? x,(t) = 


* From here on the symbol { Xq(t)} ¢ denotes the first element of the column matrix within the braces, 
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Then x(t) = Exa(t) would be a periodic function satisfying (10) and the cross-over Conditions (12). But in this 


case, the values of X(t) at times ty, tg, ... , t, might not agree with the given abscissas of the breaks or jumps of 
the characteristic X(t) = Xq given (q=1,2,...,h) or with the notation System (11) 


> given (29) 


We now determine ty, tg,... , t, starting from Requirement (29).* Going over from function x(t) to x(t) 


with the help of (7), Condition (29) is written as 


= 2, givent Pati 


{Xqo(ti, te,.- 


= Zagivert (a= 4, 2, +h). (30) 


Equations (30) serve to determine t,, tg,... , th, and are therefore the required period equations, 


Thus to find the unknown times ty, ty, . . . , tit is necessary to solve (27) for all Xqg, that is to obtain the 
dependence of each X¢g on all unknown passage times of separate portions of the characteristic, and to expect 


Conditions (30) to be satisfied. 


The operation of solving System (27) for the unknowns Xqg_is reduced to the evaluation of the matrices 


inverse to 


and Ng, and to the evaluation of their derivatives. For large h and n these computations are cum- 


bersome but the problem of construction of period equations is basically reduced to the solution of simultaneous 
linear algebraic equations and integration of the basic differential equations is not required, 


When h = 2 (two parts characteristic) it is convenient to solve Equation (27) graphically, and also when 
h = 4and the behavior is symmetrical, Then the number of unknown times is reduced by half (see [4]), the last 
Equation (27) is written in the form Xh/s,2 = ~Xqg, and the total number of equations is also reduced by half. 


We write Equations (27) showing explicitly all the unknowns 
H(t) the elements of matrix Nq(t), analogously to our notation Re 


Denoting by 
by Qf) ¢t) elements of column matrix S,(t), we obtain 


Rox (ty) + Ris (ts) + 
— Hy (+ O) — Hy, (+ 0) xi2—. 
Rox + Rin (t1) + 
— Ho (+ — Hy (+0) —... 


(t:) 202 + RY (ts) 

— Hy” (+ 0) — Hi? (+0) —.. 
Roa (te — th) + Ria (te— ts) 

— (+ O) — (+ 0) —.. 
Ri” (th — 

— Hy? (4+ O)ate+... 


Xe) for elements of M(t), and 


+ — 


O) = Qi1(+ 9), 
+ Rn-1,1 — 


+ (4) — 


0) -(+-0), (31) 


+ 2 (te — th) 


Ay-1,2(+ O) x33 Q2(+ 0), 


+ (tn — tra) 2 — 
+ (+ O) *=Qf (+0). 


The matrix of System (31) is a matrix of 2n truncated rows; it is of order hn. 


* By this requirement the q-th stage ends just as X(t) reaches the value Xq given- In a different case, for example, 
in a system with backlash, the q-th stage ends with velocity X(t) becoming zero, For cases of this kind the con- 
ditions for ending a successive stage of motion can be generalized, as shown in [1]. 
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When computers are used to calculate the solutions of (31), one can first assume arbitrary values tg, then 
rar the algebraical equation to find the unknowns ae and subsequently to look for such tq that all Xe = Xa gt gea* 
Pa +1 (q= 1,2,...,h). This cam be done automatically by introducing an aggregate parameter,say I(t, 


— Te given— and then by minimizing f(t, tg, «+ « » th). If periodic solutions exist, 


q=1 
there are also values 4, t,... , %, making the aggregate parameter equal zero. 
Having found values ty, tg,..+, hs should be verified whether they satisfy two conditions: 1) the inequal- 


ities &<...< thi should be true; ) function x_(t) constructed with the obtained ... , t, should at 


t = tq reach the value x given * P + 1 for the first time, The values of ty, ty, ..., th, which do not satisfy even 
one these conditions be re 


5. Detailed Remarks on Computation 


Computation of the derivatives of function Ryg(t) may involve difficulties due to the fact that the Fourier 


series differentiation weakens convergence; divergent series can also occur, To avoid these difficulties we trans- 
form (16): 


1 
Str 2 (s = ira). 
There > denotes the sum with the omitted one term corresponding to r = 0, 
The easy to verify identity 
Lq(s) = + gh + gist... + 
is used to obtain 
Hence 
The continuous part gives a discontinuous saw-tooth function 
after k successive differentiations. 
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The series cannot be differentiated further,term by term,as it would then become divergent; the function 
2 a . = can, however, be differentiated directly, Its derivative, equal to — =, is continuous everywhere except 
at points ¥T (integral v), where it is undetermined, The values of the derivative, however, on the left and on 


the right of point YT are equal; that is at these points the derivative is formally continuous. 


The other parts of Ryg(t) are continuous together with their derivatives up to the (n —1)th order inclusive. 
This can be deduced from the form of their Fourier coefficients (see [5]). 


We introduce functions Pr (t), 
(33) 
and (t), f(t), (t), 
(34) 
defined as follows: 
1 
P2(t) = (ire)? 


t 
i 
iy 


1 eirat 


O<t<T) 
and so on ;* 
4 i n—2 


f(t) dt +a 


4 


rT 
* It should be mentioned that for integral k > 0 


— 1) 


(iro 


where ¢(k) is Riemann zeta-function tabulated in [6], page 444. 


L, (ire) 
n— : 
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t 


and so 


Having now written all functions RE) (x) (j= 0,1,..., m1) we note the fact that everyone of them can be 
represented as a linear combination of Functions (33) and (34) with coefficients gt For 
example 


(t) = — (t) — (t) —. — 


This fact simplies considerably the computation of the values Rf) (r). 


Coming now to the problem of computation of Ng (+ O) and Sq (+ O), we agte that due to (28) the question is 
reduced to finding a method of computation of M, (+ O). When the elements RB of matrix M(t) are constructed 
by means of Functions (33) and (34) it is obvious that substitution t = 0 gives a correct value of matrix M_, (+ 0). 
One should especially be careful not to let argument t of matrix M(t) leave the interval (O, T) as Functions (33) 
are defined only in that interval. 


6. Possible Generalization 


Finally, we should note that by using the method given in this paper one is able to reconstruct the periodic 
solutions of a wider class of simultaneous differential equations, which have on their right-hand side periodic 
functions F;(t), period T = 


= + + Fil), 
j=1 


y=f(%,) (=4, 2,...,A). 


When the functions F,(t) are smooth and differentiable the required number of times, the derived Equation 
(2) is of the form (see [1]): 


D(p)z = K (p)f (2) + O(t). 
Equation (10) becomes 
Lq(p)z(t)=D(t) (tg 
Simultaneous Equations (26) and (27) for finding the unknowns X(t) and Xgg becomes as follows: 
Xq(t) = Mg (t — Xq—1,2 — (t — te) — Sq (t — tq) — Py (t), 


M, (t, — to) — Ni (+ O) Xin = Sy (+ O) + Py (ty), 
My (tz — ty) Na (+ O) Xog = Sg (+O) + Po (ts), (36) 


(35) 


= Xe: 


On right-hand sides we find now additionally column matrices 


* We note that 


(irw)"—* (irw)"—* 
L, (ire) ~ 2 (ire) 
as 


L,(—ire) “Ly (ire) * 


Mx (th — ths) — Nn (+0) Xna = Sn (+0) + Prite), 
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T, (t) 


(t) 


Te» (t) 
Function T ,(t) is given by the formula 


ta 
r 


As distinct from (26) and (27) the quantity w= 2n/t, is known in (35) and (36); there appears, however, a 
new unknown t, which in the self-oscillatory case was equal zero as we were then free to choose the beginning 
of the time scale. 


APPENDIX I 


Derivation of Crossing-Over Conditions 
In agreement with (7) when 
Zz (t*) = (t*) + Pq 
y = kee (t*) — = (t*) + + 


and when 


(t**) = 2 (t**) + 
We form the differences 


(t**) — = (t**).— (t*) + — Pg 


y (t**) —y = ke + + Rott ae Aq 
and then we make t®® and t* approach tq from the right and left,respectively. In the (11) notation we obtain 


bo = — + Pata — Par 
Applying the same limit procedure to the differences of functions x(t) and y(t) we determine € , and ny, (k = 
=1,2,...,m-1). Substituting those in (3) we obtain after some rearranging 


Solving these simultaneous equations for all x‘,(k = 0, 1,..., n—1) we shall obtain exactly System (12), 
where Ugk and Vak are determined by the recurrence Formulas (13), 


APPENDIX II 
Construction of Functions Rygl) 


We introduce generalized differential operator p* (see [4]) which is defined in the following way: 


p*F (t) = pF (t) + >) (t— ty), 
@ 


(t) = (t)+ >) + tes 
q 
(t) = (t) + (t — t,) + +t) 
q 


etc. 


There p = d/dt is the ordinary differential operator; ty ry Effetc. are the magnitudes of gaps of the piece- 
wise continuous function F(t) and its derivatives dF/dt, @F/de etc, at time t = tgi 5(t), 5 (t) etc, are Dirac impulse 
function and its derivatives, | 


For the function Ryq(t) according to Condition 2 in Section 3 we have 


Rig (t) = (t)s 
P*R (t) = PRiq (t), 


In agreement with Condition 3 (Section 3) 


By Condition 4 (Section 3) 


() = + 8 (t—vT), 


(t) = (t) + —vT). 


v=——oo 


Multiplying (II. 1), (II. 2) and (II. 3) by gf, gd ayrtete gs, the coefficients of operator L,(P). and adding 
on the left and right-hand sides we obtain 


Lg (P*) Rig (t) = Lg (P) Reg) + 4-1 


%) 
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If 


(w= 2) 


then, using a known formula in Fourier series theory, 


wo 
and also from known formulas on Fourier analysis of periodic impulse function) 8(t—vT) and its deriva- 


tives 
@ 
8” (t — vt) = > (irw)%e'™ etc. 


and by taking into account Condition | we obtain 


L, (iro) y, = + +... + 


( = ire). 
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SUMMARY 


The paper deals with determination of periodic behavior of control systems having an arbitrary piece-linear 
characteristic, The said periodic behavior is determined as a complete Fourier series without neglecting harmonics, 
The problem is reduced to solving simultaneous transcendental equations that determine motion time in each part 
of the nonlinear characteristic. Some problems of using computers to solve period equations are considered. 
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CONCERNING THE EQUIVALENCE OF PULSE AND CONTINUOUS=-DATA 
CONTROL SYSTEMS 


V. A. Rubtsov 


Abstract 


The insufficiency of the generally used equivalence criterion of pulse 
and continuous-data control systems is established. Sufficient conditions for 
the equivalence of closed loop control systems are determined. 


It is usually assumed that the analysis of a pulse control system can be reduced to the analysis of the cor- 
responding continuous-data system, if the repetition period T, (equal to the time interval between two consecutive 
error measurements) is small compared to the basic (essential) time constants of the linear part of the pulse con- 
trol system [1-4]. However, it will be subsequently shown that considerable differences can exist between the closed 


loop, pulse and continuous control systems under comparison, both with respect to transient processes as well as stabi- 
lity conditions, 


Derivation of Sufficient Conditions for the Equivalence of Linear, Closed Loop, 
Pulse and Continuous-data Control Systems 


Let us examine a series of preliminary proposals which will permit us to determine the sufficient conditions 
for the equivalence of the pulse and continuous control systems being compared, 


Fig. 1. 
1, Sufficient conditions for the equivalence of closed loop, pulse and continuous control systems cannot be 
established from the comparison of the transfer functions of the open loop control systems, 


A linear, pulse control system usually includes its own linear part, with transfer function K(p), and a pulse 
feedback loop with the amplification factor of the pulse unit Kp and the parameter 7. 


Blocks schematics of the closed, and open loop, pulse control system are shown in Fig. 1, where the follow- 
ing designations are introduced: Ke (q) — transfer function of the closed loop control system, K*(q) — transfer 
function of the open loop control system, y — a parameter characterizing the mark-space ratio of the error 


measurements (of the succession of pulses), 0< y= 1, Kip) = GES + Mp) and Q(p) are polynomials in p. T, will 
be used to designate the pulse repetition period, 


Let us express the transfer function of the linear part K(p) in the form of a fraction-rational function, making 
at the same time, the substitution p = q/T;: 


peo Kig Your 
Y 
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Bog™ +B. (1) 

K (q) = ke mf _+ + + bra + bo 
+4,_,9" 419 + dy 

Here a, = 1, a; are functions of the systems time constants T, and repetition period Ty, ky 1s a constant coef- 

ficient (the amplification factor of the linear part). 


Expression (1) can also be expressed in the form 


n P( 


where qy = T,py» Py — are simple poles of the transfer function of the linear part of the control system with a pulse 
feedback loop. 


If the condition 


lim K (q,) = 0, (2) 


holds, i.e., in Expression (1) the degree of the numerator m is less than the degree of the denominator n by at 
least one, then the transfer function of the open loop control system, including the pulse unit, can be written in 
the form 


= P i— et (3) 
(9) kokp 1 Q’(9,) qv ef 


The mathematical expression for a practical criterion which makes it possible to consider a pulse control 
system as an ordinary system responding to a continuously measured error, is usually written in the form [2] 


Tn Pv\ = | @v| << 4. (4) 


The physical significance of Inequality (4), as it is well known, is that all, or at.least the basic (essential) 
time constants of the linear part of a control system are large as compared with the repetition period T; [2]. 


Moreover, Expression (4) emerges as the result of the establishment of equivalence between the transfer 
function of an open loop, continuous-data control system (1) and the transfer function of an open loop pulse 
control system (3), for the limiting transition, i.e., for Ty; > 0. 


In fact, Condition (4) is necessary and sufficient for the establishment of equivalence between pulse cir- 
cuits and corresponding continuous-data circuits, if one considers the usual pulse circuits which do not contain 
a feedback loop, Because the poles of the transfer function of the open loop, pulse control system coincide with 
the poles of the transfer function of the linear part of the system, or differ from them by #2 kj, i.e., 


then in the region s < Im qys # the minor poles qj will always correspond to the minor poles qy (when 
lavl << 14. ¢ y, for T; — 9, when yl<< 1), it is permissible to make, in Expression (3), not only 
the substitution eY = 1+ qy, eW7= 1+ ydy, but also e4 = 1+ q, It is not hard to see that, after this substitution 
in Expression (3), the transformer function of the open loop pulse control system will differ from the transfer 
function of the linear part of this system only by a constant multiplier k.y. However, the insufficiency of Cri- 
terion (4) for closed loop pulse control systems becomes obvious if one a few concrete examples, 


Let us examine the simplest static control system using a nonperiodic element with a time constant T. 


tive 
closed 
} stabi- 
o v= Wa (5) 
king, 
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Let K(q) = ky as 


» Where 8=1,/T and y =1. The behavior of the closed loop control system is deter- 


mined by its transfer function and, partially, by the poles of this function, Making use of (3), the pole of the 
transfer function of the closed loop control system can be expressed in the form 


= + In [e-* — A(t —e-4)], (6) 


where A = kgkpy- Correspondingly, the pole of the transfer function of the continuous-data system will be 


Obviously, the equivalence of the systemscan only be considered when 


I: (8) 


Let us suppose that Condition (4) is satisfied, in which case (6) can be written in the form: 


qo In[{1 —B(1 + A)]. (6") 


If we follow the commonly accepted approach [2-4], then Condition (8) must always be satisfied if 
Inequality (4) holds, However, as can readily be seen from (6") and (7), this will only be true when Ag << 1, L.e., 
for a sufficiently small feedback factor A (small amplification in the open loop pulse control system). 


From (6°) and (7) it follows immediately that, for 2> Ag = 1, the systems differ considerably, even when 
Condition (4) is satisfied, Thus, for example, in the continuous-data system, the transient process will be non- 
periodic for any limiting values of Ag, while in the pulse control system, for 2> Ag = 1, the transient process 
will be oscillatory. Moreover, the pulse control system will already be unable to control when Ag * 2, because it 


will be self-exciting. 


I 
| 
bat cal 
In 


42 00 100 A*Ky hy 
Fig. 2. 


What has been said is illustrated in Fig. 2, which shows a curve 
of the function n = f(A), when n is the duration of the transient process, 
expressed in repetition periods . oa I is the curve showing the pehavior 
of the pulse control system when g = 0.05 andy = 1, II is the curve 
showing the behavior of the compared® continuous-data system, The 
region in which the transient process in the pulse control system is of 
an oscillatory nature is shown shaded in Fig. 2, In determining the 
duration of the transient process n it was assumed that the transient 


process ends when gy) —Uout = gor Where gp is the 
stable value, 


Thus, from the examination of the above example of the 
simplest form of control system it becomes apparent that Condition 
(4), which results from the analysis of compared open loop, pulse 
and continuous~data control systems, is not sufficient for the equiv 
alence of these systems in the closed loop condition. Inequality 


(4) is, in the general case, a necessary but not sufficient condition for the equivalence of closed loop automatic 


control systems, 


In the examination of closed loop, pulse control systems, the behavior of the system is characterized by 
the transfer function of the closed loop system, which can be written in the form: 


= 


(q) (9) 
1+ K (q) 


* By a compared continuous-data system is understood such a continuous control system into which a pulse control 
system degenerates for the limiting value of the feedback factor, when T, —> 0. 


where K®(q) is the transfer function of the open loop system, determined by Expression (3). 
For the compared continuous control system, the transfer function can be expressed in the form: 


Kel) (10) 


Obviously, the compared closed loop systems will be equivalent under the condition that the denominator 
of Expression (9) will degenerate into the denominator of Expression (10). From the formal identity of expressions 
K*(q) and K*,(q) as well as K(q) and K,(q),{t is obvious that the conditions for the equivalence of the closed loop, 
pulse and continuous-data control systems will be analogous to the conditions for the equivalence of the open loop 
systems, 


In establishing the equivalence of the open loop, pulse and continuous-data control systems it was assumed 
that it was necessary to fulfill Conditions (4) and (5). 


Therefore, for the equivalence of the closed loop, control systems, two analogous conditions must be satis- 
fied, namely: 


Gr = Sy | (11) 


Conditions (11) are the determining ones in the establishment of a practical criterion for the equivalence 
of closed loop, pulse and continuous-data control systems, 


In fulfilling (4), the denominator in (9) can be expressed in the form 


= @’(a,) — 9) + 


n n (12) 
and, correspondingly, the denominator of (10) can be written in the form 
v c = (9y) % — qv v 


then, obviously, it is permissible to make, in 
Expression (12), the substitution® e% = 1+ After this substitution, expressions ) and G.(q,) coincide 
completely, and the behavior of the pulse control system will not differ from the behavior of the continuous-data 
system; in other words, the systems will be equivalent. 


In this way, Conditions (11) will be necessary and sufficient conditions for the equivalence of closed loop, 
pulse and continuous-data control systems, However, Conditions (11) are of little practical value since they 
require the determination of the numerical values of the poles of the closed loop, pulse control system's transfer 
function, 


In practice, it is desirable that the sufficient conditions for the equivalence of pulse and continuous-data 
control systems be completely determined by the characteristics of the continuous control system itself, and do 
not require the direct analysis of the pulse control system, It should be noted that, for closed loop control 
systems, Criterion (4) would be sufficient only in the case where, from the condition | q,|<< 1, condition lag yl<< 1 
would always follow. However, such a mutual relationship does not occur in general. 


*In paper [i], in establishing the equivalence of pulse and continuous control systems, it was postulated that the 
expansion ede = 1+ q. was permissible, However, such a substitution is far from always possible. 
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2. The minor roots of the equation G_(q_) = 0 correspond to the minor roots of equation G* (q*.) = 0. 


In the case where the linear portion of the pulse control system can be expressed as a fraction — rational 

function of the form (1), whenever Inequality (4) is satisfied, there will always be a definite relationship between 
the roots of the equations G.(q,) = 0 and G* (q%) = 0, i.e., between the poles of the transfer functions of the 
closed loop pulse and continuous-data control systems, This relationship can be expressed in the form 


+ q,). (14) 


Expression (14) results from the fact that G&(q%,) can be formed directly from G,(qc), whenever (4) is satis- 
fied, by means of the substitution q, = z —1 [2] in G.(dc), where z = et, Expression (12) was obtained exactly in 
this way. 


If, in (14}, | qgy] << 1 then |q@y | will also be quite small and the approximate equality of the poles 
results, which can be written in the form 


Gey = (1 EF) (15) 
where 6 << 1, 


From what has been said above it follows that it is sufficient to determine the condition for the transfer 
function of the closed-loop, continuous~data system to have small poles, to thereby establish the required suf- 
ficient conditions for the equivalence of the compared pulse and continuous-data control systems. 


3. For certain values of the coefficients in the denominator of Expression (1), the roots of the denominator 
will have small values, i.e., |qy|<< 1. 


The coefficients of the denominator of transfer function (1) are related to its roots thru Viet's equations, 
This relationship can be expressed by means of the following equalities: 


a, = i, 
= — (91 + 92 +--+ -+49,)) 


(16) 
On—g = — (919293 + 919294 + 919994 In? 
Gy = (—1)"9192- - - 
From (16) the following inequalities follow: 
An—1 < n nex 
< (n — 1) (9,93) max 
(17) 


— 2) (9,995) max 


where q; 19, 18 the lar,est root, (q44j)max 18 the largest product of a pair of roots, (4;4j9; max 18 the largest 
product of three roots, etc, We will conditionally assume those roots to be small for which the following relation- 


ship holds 
| q, | < (18) 


<9? max 
930 


where o << 1, 
If (18) is fulfilled then Inequalities (17) can be replaced by the following conditions: 


Ge, 

(n — 1) 0°, 

(29) 
ay <0" 


From (19) it follows that if the roots of the denominator of Expression(1) are small, its constant coefficients 
are also small or, at least, cannot be greater than the values determined by Inequalities (19). 


Let us limit the class of polynomials expressing the denominator of transfer function (1), to those polynomials 
which satisfy the Hurwitz conditions, In that case it turns out that, if Conditions (19) are satisfied, the roots of the 
denominator of (1) will always be small compared to unity and on the other hand, if the values of the coefficients 
do not satisfy these inequalities, then the conditions for smallness of the roots are not fulfilled, at least for one or 
several roots, 


4, Sufficient conditions for the equivalence of closed loop, pulse and continuous-data control systems are 
determined not only by the systems time constants T,, but also by the feedback factors A. 


Taking (10) into account, let us write the denominator of the transfer function of the continuous-data system 
being compared in the form 


Geldc) 1+ (20) 
In accordance with (1), for a, = 1, the right side of (20) can be expressed in the form 


(dna + (dy + Abs) Ge + + Ady. 


It follows from (20') that the coefficients of the denominator of the transfer function of a closed loop control 
system are functions of not only the time constants of the different loops in the system and the repetition period T, 
but also of the parameter A. In present pulse control systems, the value of A can reach thousands and even tens of 


thousands and, at the same time, the coefficients of Equation (20°) can not be small, even for sufficiently large 
time constants of the system, when Criterion (4) is fulfillable. 


In the investigation of a closed loop control system, Conditions (19) will be satisfied, in general, only for 
certain values of parameter A. In accordance with (19), these values of A will be determined by the following 
inequalities; 


(21) 


where o << 1, 
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The roots of the denominator of the transfer function of the closed loop, continuous-data system (10) will be 
very small compared to unity specifically for those values of A which satisfy Inequalities (21). Consequently, to 
the extent that (15) will be satisfied at the same time, Conditions (21) will be the sufficient conditions for the 
equivalence of the closed loop, pulse and continuous data systems, 


The degree of the inequality 0 <<1 determines the accuracy of the coincidence of the basic characteristics 
of the transient processes and the stability conditions, and the greater this inequality, the more are the compared 
systems equivalent, For most practical cases it is sufficient to take o = 0.1, Thus, for example, in the example 
investigated above, when o = 0.1, the difference between the durations of the transient process in the compared 
systems will not exceed 1%, 


It should be noted that complete equivalence of the systems, i.e., their equivalence not only with respect to 
the nature of the transient processes, but also with respect to stability, requires that the sufficient conditions be ful- 
filled even for the value of A corresponding to the excitation condition, i.e., for A = Aq. However, the compared 
systems can be equivalent with respect to transient processes by fulfilling sufficient conditions not depending on 
the fulfillment of the conditions for total equivalence. An example of such a system is a continuous-data control 
system, described by a third order differential equation, for which, in Expression (1), the coefficients b,_y = b,-,= 0 
and by = a. 

In fact, if the linear portion of a control system, having a pulse feedback loop, represents a three-section 
filter with equal time constants, then the feedback factor A (amplification factor) for which the system becomes 
excited will be equal to 8 [5] irrespective of the absolute value of the time constants. Hence, by increasing the 
system's time constants one can always satisfy sufficient Conditions (21) for A = Ag, i.e., to have the compared 
systems completely equivalent, 


However, if one of the coefficients b,; 4 or b, __g is not equal to zero, then the continuous-data system 
being compared can, as it is well known, remain stable for any large value of A [5]. Consequently, in this case, 
the condition for complete equivalence will not hold and the conditions for the stability of the pulse control system 
can be determined only by the direct analysis of such a system. 


SUMMARY 
1, The condition for the equivalence of pulse and continuous-data control systems 
is necessary but not sufficient in the investigation of closed loop control systems, 


2. To establish the equivalence of the systems with respect to transient processes it is necessary that, besides 
Condition (4), sufficient Conditions (21) be satisfied. 


3. The investigation of the stability of closed loop, pulse control systems, based on the analysis of the cor- 
responding continuous-data systems is, in general, not possible. Such an investigation is possible for a particular 
class of closed-loop systems for which sufficient Conditions (21) are satisfied even for A = Ag. 
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CONCERNING THE EXISTENCE OF A CYCLE BEYOND THE ABSOLUTE 
st to STABILITY CONDITIONS OF A THREE DIMENSIONAL SYSTEM®* 


pared B. V. Shirokorad 


9 Abstract 


n Beyond the absolute stability conditions (according to Lur’e-Letov [1]) of a 
es closed-loop system with a neutral object and a regulator having a nonlinear speed 
ne of control reset, there always exists a nontrivial (in particular, periodic) stable 

od regime. Physical interpretations of the system are presented, 


A study is made of the properties of the phase portrait of the simplest type of automatic control system [2] 
described by the equations 


—% +f (2), 
= + Bota (9), 


for the critical case of a neutral object p; = 0 when 6; = Bg > pg > 0. 


If, in a finite interval, the function f (0) belongs to A** class, if it is aLipshits-type and is nonconvex 
(part 2), it satisfies the conditions of dissipativity (20), repulsivity (26), and Bellman's condition (27)* * * , then 
there exists a sphere S and a parallelipiped P situated inside it which contains a closed-loop trajectory-cycle 
(stable regime, part 6), Trajectories which get into P at t = 0 are contained in S when t> 0 (are stable according 
‘or- to Lagrange) and return into P whent> T> 0 (mainly dissipative, according to Massera, part 2). Any trajectory 
lar passing thru a point in P and not tending toward the origin as t—> + a, has among its w-defined set [3, 4] an 
oscillatory regime (part 4). In the vicinity of the origin there exist only two unstable, in the Liapunov sense, 
trajectories (part 3), tending to the origin as t—> + o. 


In the analysis of automatic control systems, in particular (1), the basic problem is the determination of the 
steady state, stable regimes(part 6) and the investigation of transient processes characterizing the method by which 
other (transient) regimes approach these stable ones, 


In the conditions for total asymptotic stability By < p, fiyax > a, when |o| —> o [2], the phase 


portrait of System (1), in the case of a neutral object, does not contain any closed trajectories (nontrivial, stable 
modes), All the other (transient) regimes, which fill the portrait, approach the origin (zero solution or trivial, 
stable regime). For Bg > pg there always exists in the A class such an f(0), for which the portrait contains at 
least one trajectory which leaves a sufficiently small region around the origin for a sufficiently large t > 0 (part 


[957 


and * Delivered on October 9, 1957 at 1.A.T., AN SSSR during the seminar on automatic control, 

**To the A class belong certain, piecewise-continuous, with a finite number of discontinuities of the ist type, 
apriori fixed functions f(o), satisfying the condition f(0) > 0 when | o] > o, = O and f(c) =0 when fo |x Oe. 
For the sake of simplicity, henceforth 6, = 0. 

***For example, if the function f(0) is piecewise-linear, 
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3), At the same time, it might turn out apriori that the portrait contains some finite region, containing the origin, 
from which the inside trajectories do not emerge for t > 0, and into which enter all the outside ones for sufficiently 
large t> 0. If the diameter of this region is so small that the technical requirements on the accuracy of control 
are satisfied, then recommendations concerning the adjustment of the control unit, based on the theory of absolute 
stability, can become incorrect. 


These and many other factors make it interesting to study not only the local properties of the zero solution of 
(1), but also the properties of the phase portrait as a whole. 


The framework of the problems of B, V. Bulgakov [5-7] is part of the work started back in the ‘40's by A. A, 
Andronov [8, 9], It complements the work of A, I. Lur'e [1] and A, M, Letov [2, 10] and, on the other hand, 
generalizes the work of Giuseppe Colombo [11]® for the case when the coefficient of coupling between the grid and 
plate circuits a vanishes (see § 2, (6), [11]), which corresponds, in our case, toCondition (7) (neutrality of the 
object py = 0, pg > 0, see further on), while the function f(0) belongs to the A class, where atte its smoothness 

nor its symmetry are presumed, 


In System (1) being studied belongs to a general class of controlled systems, examined in [13], to which, 
in particular, belong the systems investigated in [14]. The simple apparatus used for the study of the phase por- 
trait of a three-dimensional system in [11] and in this work, can be applied to other types of equations, to those 
studied for example in [15], or to systems containing several nonlinear functions, 


Probably the statement of the problem dealt with in this work, concerning the investigation of the phase 
portrait of (1) beyond the necessary and sufficient (if one excludes from the investigation the bifurcational case 
of the boundary) conditions for absolute stability [16]:* * 


Pi = —(PiP2 + Po=—(Bit+Bs), P =pips, =p, + pr. 
is original, 


1. Physical Interpretations 


Let us consider in our analytical investigation two physical interpretations of System(1). The first is an auto~ 
matic control system [2], consisting of an aeroplane (Fig. 1) and a stabilizer (autopilot) controlling the bank angle 
¥ (rotation about its axis of symmetry OX,). The second is the classical circuit of a single tube oscillator (Fig. 2) 


[11]. 


Turning to Fig. 1, let us consider the simplest system for the automatic piloting of a'plane, flying along a 
straight course, at a constant velocity, at a given height above the surface of the earth. We will assume that the 
channel for stabilization with respect to the OX, axis depends so little on the two other channels (with respect to 
the OY, and OZ, axes) that the absolute angular velocity wy = y of the ‘plane's rotation about axis OX, depends 
only on the deflection 6 , of the aileron (bank control) from its neutral position. 


*The basic theoretical principles of [11] and the present article can be found in the popular report of S, A. 
Stebakov [12]. The author of [11] admitted a series of inaccuracies into his paper and even an incorrect evalua- 
tion of the value of y (see (24) and above, in paper [11]. 

* *in a two-dimensional, speed and position control plane, the quantities py and pg take on a very definite physical 
sense: py, is the position transfer ratio of the stabilizer, p, is the velocity transfer ratio. 


The proof of these conditions [16] required that P > 0; however, this limitation is not essential since, in our 
case, whenp, = 0, 0, they coincide with the previously known gy, < 0, Bg < pz [2], which guarantee sufficiencye 
The necessity is proven by the function of the A class f (@) = Do, D = const, satisfying the repulsivity Condition (26) 
and upsetting the stability when B, > ps. 


le 
2) 


Fig. 1. 


Introducing into the analysis, in our case, the constant, 
negative, aerodynamic coefficients My* kg sec® of natural 
damping, and Mea kg of the aileron's effectiveness, we can 
describe the banking motion by the equation 


where I, kg sec* is the constant moment of inertia of the aero- 
plane about its principal axis OX;. The bank autopilot, with 
strong feedback coupling, is described by the equations 


where the positive, constant quantity ky is called the autopilot's transfer ratio, while the class A function f(o) is 
the speed characteristic of the steering machine, 


The system of Equations (2) and (3) are reduced to the form of (1) by the following substitution: 


Fig. 2. 


— MSs 
4 
(4) 


where T is a dimensionless time parameter, 


Turning now to the circuit of the single tube oscillator (Fig. 2), we will use i, V, and V, to designate the 
triode's plate current, plate and grid potentials, respectively. If p is the amplification factor, 1= 9(Vp + 
+ Vg), where the continuous function ¢ (of the A class) is called the triode characteristic (Fig. 3). 


Furthermore, let E and Eg represent, respectively, the sources of emf in the plate and grid circuits, C,2, and 
1 the capacitance, inductance and resistance in the grid circuit; L and R the inductance and resistance in the plate 
Circuit; iy and ig the currents in the grid circuit, passing thru the resistance r and inductance! . 


If M is the mutual inductance between the plate and grid circuits (feedback factor) then the oscillator equa- 
tions take on the following form: 


E, —li,; — Mi, 


= i, + (5) 
V, = E— Mi,—Li— Ri, 
=Vp+uVeg, 


in, 
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te 
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where g(i) is an essentially smoothly growing function, namely, the triode's reverse 


is characteristic (Fig. 3). 
Let us examine the system of equations in (5) for the particular case when the 
inequalities 
<= 
Fig. 3. (6) 
and the equality 
(7) 
M 


are satisfied, which is always possible by proper selection of parameters, 


To reduce (5) to the normal Cauchy form of three equations, we will introduce in Conditions (6) and (7), 


the designations 


(8) 


and examine the equation ini 


Ri + g (i) =E—pEy, 


which has the single solution i = 4, > 0. This solution determines the steady state value (stable or unstable) of 
the triode’s plate current. 


Let us further designate the deviation of the plate current i from the quiescent iy by means of 


and introduce into the analysis the function 


= + (io + 2) — 8 (10) 


which vanishes together with the deviation x and grows smoothly, since it follows from (6) that 


If we now introduce into the analysis the new variables 
y=z+/(2), (11) 


then, taking (8-11) into account, System( 5) can be reduced to the normal form 


(12) 
Zz == — — k,f (2), 


to which System (1) also reduces by means of the substitution 


y= + Boe, 
z= —k,2, (13) 


(6) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


—B, = = — 


pi = 0, > 0, Ps 


2. Lagrange Stability 


For a Lipshits-type,* * nonconvex® * * f(x) of the A class, under Conditions (14), is examined in Systems (1), 


expressed in the form (12) belonging to a class of dynamic equations [17], because the right sides of its equations 
do not contain explicitly the time parameter t. 


For the sake of brevity, we will designate the integral curve (trajectory) passing thru the point X(0) ={ x(0), 


y(0), 2(0)} att =0, by X(t) = X[t, X(0)], where X(t) = { x(t), y(t), z(t} , (X(0) = XL[0, X(0)]) is a vector describing 
the point's motion, 


If there exists a sphere S, with radius R, within which is contained the whole trajectory X(t), i.e., | X(t)] < 
< R(| X(t)] = GF + y+ 27) for every t, then the trajectory X(t) is called stable according to Lagrange; if this 
condition only holds for t> 0(t< 0), then X(t) is called stable according to Lagrange in the positive (negative) 
direction, 


Let us prove that there exists such an open parallellipiped P 
(15) 


(where € and € are some quantities) that every trajectory X(t) of System(12) which passes through the point X(0)€P 
at t = 0, does not emerge from the sphere S(X(t) €.S) when t > 0, L.e., it is stable according to Lagrange in the 
positive sense, Moreover, there exists a number T > 0 such that X(t) € P whent> T, or System (12) is mainly 
dissipative (belongs to Massera‘s D-class [18]). At the same time each trajectory X(t) = X[t, X(0)] (X(0) © Pjis 
determined for an infinite time interval, since afortiori the following theorem [3] takes effect: if the trajectory 

of motion remains, with increase in time, in a closed, bounded region, then the motion can be extended over an 


. infinite interval [t), + 


For the proof let us examine the system 
V+ py thy tay—/(2), (16) 


which is obtained from (12) if z is excluded, 
Let us suppose apriori, that the function x = x(t) is known and the conditions of absolute stability will be 
k 
violated, i.e., — 
> Ps 


Then 


y= + C,sin t)exp(— ) + 
+ (2) sin — t) exp — t)dt, 
0 
* Depending on the value of the mutual inductance M, the considered oscillator represents, as it is shown [see (13), 
(14)], not only an electronic model of the 1st circuit (with the restriction (7)), but, as it can easily be seen, the 


2nd circuit in B. V. Bulgakov's problem [2]. Here, M< L/y corresponds to a stable object with positive self- 
recovery, while M > L/y corresponds to an unstable (with negative self-recovery) one. 


**For a sufficiently small | h| :|f (x + h)—f(x)| < const h, These Lipshits conditions can be replaced by Osgoode's 


conditions, or other less restricting ones, guaranteeing single-valuedness and continuous dependence on the initial 
conditions, 


***F| [Oxy + (1—@) x J] < (xp! + F when 0<6 <1 for all | x| and | . 


only if* 

| 

| 

- | 
| 

= 
| 

| 
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= — pt = Bky + — px) > 0 


C,=y(0), = (0) + 22 (0). 
Integrating by parts, we obtain 


y ={y(0)cos-5-t ++ {pay (0) + 22 (0) — (0)] sin t} + 


t 
+ {[cos (t sin 3 (¢—*)oxp[ —9)]} dt. 
0 


Differentiating both parts and noting that ¥ = z, we obtain 


= {[2(0) — (0)} cos + (0) — 2kyy (0) — paz (0)] x 
a 
x sin > i} exp(— + — sin (t—1) + 


+ pcos (t — (x —t)}} dt. 


Let us suppose that the starting point X(0) belongs to a sphere C, with center at the origin and radius r, i.e., 
X(0)€C, and let us examine, in the time interval [0,7] the restricting inequalities resulting from (17) and (18) 


ly(t)| <<mr+kz, (19) 


k= 2k, (> + 


+ 2k, + es 


n=1+kh,+ bo k,(34+ Z—maxjz(t) 


Let us suppose that the class A function f(x) satisfies the following “dissipativity” conditions. There exists 
a sufficiently large r > 0 such that the roots of the equations 


ma =f no+ly=/ (20) 


corresponding tow=r, and € ,9<0when w=~-t, lie in the interval r] and moreover, the following 
inequalities are satisfied:* 


(21) 
* Inequalities (21) are satisfied, for example, for the function f (€), for which + when —> +o. 


In fact, applying "Hospital's rule, we obtain 
lim / r = lim m/ |f’ — n] = Owhen r + 


$ 
where 
and 
(18) 
where 
4 
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Then the following inequalities are fulfilled: 


x(t) <§ whent >0, (22) 


only if E< x(0)< €. If this is not fulfilled then two cases can arise, In the first there isa minimum Tf > 0, for 
which x(T) = € and X(T)2= 0, However, in this case, a contradiction results because 
x(t) = y(t) — < mr + =O. 
A similar contradiction arises in the second case when there exists a minimum T > 0 for which x(T) =§ 
and X(T) < 0. 


From Inequalities (19) and (22) it follows that every trajectory of (16), passing at time t = 0, through the 
point X(0), lying within parallellipiped (15), is Lagrange~stable in the positive sense, since every such positive 
semi-trajectory (t > 0) is within the sphere S with radius 


Let us now prove the existence of a quantity T > 0 for which the trajectory X[t, X(0)] €P whent> T, if 
X(0) € P. 


There exists a number € > 0 which is such that the inequalities 


r—e 


are satisfied, which is a consequence of (21). 
Let us denote T = max (Ty, Tg), where 


From a more accurate evaluation of expressions y and z than in (19), there result from (17) and (18) the fol- 
lowing restricting inequalities: 


— mr oxp(— t) + <y ()< + mreoxp(—Ht), 
—nr exp(— 4. t) +- In<z(th<ly + nr exp(— t) 
From Inequalities (22) and (23) it follows that 

*For Lagrange-stability it is necessary that only Condition (22) be fulfilled, which is guaranteed, for example, only 
by the first Inequality in (21). All the other formulae containing z [(18) and further] are, in this case, superfluous, 
Indeed, from (16) it follows that =—p,¥ —kyf (x). From this | pp | ¥ | + kyl f(x)< pay (py + ky f ), where 


=f max f (x). According to Tonelli's lemma [19] z = § is bounded together with y. 
Es 


| 
| 
|| 
) 
0) i 
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3. Liapunov-Instability of the Zero Solution 


Certain local properties of the zero solution of System (12) will be needed later on, in particular, its repulsive- 
ness (repulsion of the trajectory in certain special neighborhoods of the zero solution). 


For the sake of simplicity, let us restrict our examination to those class A functions which, in the neighborhood 
of zero, are represented in the form 


= Dz + 9(2), (25) 
where the constant D satisfies the repulsiveness condition 


0< (26) 
while the function (x) satisfies Bellman's condition [20)” 


p(x) = o(2). (27) 
In this case, the linear approximation of (12) is the system 


= y — Dz, 
ves, (28) 
z= — pz — k, Dz, 
the characteristic equation of which 
+ (D + + Dd + kyD = 0 (29) 


has an unstable zero solution, because its Hurwitz determinant is 


H = — p2D (B2— p2 — D) <0. 
On the other hand, when the repulsivity Condition (26) is satisfied H> 0, and the zero solution of (27) is stable, 


Fig. 4. Fig. 5. 


It is known (direct and reverse theorems in Liapunov's second method) that, under Bellman's Condition (27), 
the local stability and instability properties of the zero solution of the linear approximation (28) are carried over 
to the zero solution of the nonlinear System (12). 


* o(x) denotes a quantity of an order of smallness greater than the order of smallness of x. 


940 


ive- 


5) 


6) 


7) 


8) 
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stable, 


27), 
er 


Under Conditions (26), only two cases are possible for the roots of the characteristic equation of (28), to the 
first of which correspond the three real roots A, », ¥, while to the second, one real and a complex conjugate pair 
A, w+1iv¥, wiv. The values A, yp, Y are, in both cases, positive. 


In the first case, the origin is an unstable, singular (particular) point of the “node-saddle" type (Fig. 4). 
System (12) is reduced to the form 


ey, y = vy ey, 
where €4, €g and €s are an order of magnitude smaller than (x* + y¥* + oy, 


From here, in a sufficiently small vicinity of the origin there is an open cylinder Cy, center at the origin, 
radius pp) = (2° + yyA and height E"E (Fig. 4), which is such that all the trajectories emerge from its side surfaces, 


since + xt + + 2 + yeq)> 0, and enter thru its bases F and Fy (faces), since there 
=~2\2" + Qzey< 0. The cylinder’s vertical axis EE lies along the straight line 


y=(D—)d)z, z=—by, (30) 


where D~A =— pg—y —¥ < 0, passing through the first quadrant (x> 0, y< 0,2> 0) when x> 0,and through the 
second (x < 0, y> 0,2< 0) when x< 0, From this it follows that the diameter of the region containing C, can be 
made so small that the base F will lie in the fitst, and Fy in the second quadrants, while the projections of these bases 
F, = 9 and Fyz = » onto the planes z = 0 will lie, correspondingly, in the fourth and second quadrants of the plane 
z=0, ‘ 


In the second case of the location of the roots (—A, yw + i¥, y — iv) the picture is analogous to the first case 
(Fig. 5). The origin is an unstable, singular point of the “focus-saddle" type. In a sufficiently small vicinity of 
the origin there is a cylinder Cy with properties similar to those of the first, and the bases of which F and F, are 
also situated, respectively, in the first and second quadrants, while the projections Fy = 9 and Fyy = » are in the 
fourth and second quadrants of the plane z = 0, since, in this case too, D— A =— pg ~2y < 0, 


Cylinder C, will play an important part later on, 


4. Separately Generated Surfaces of the Phase Portrait 


In phase space, let us examine two cylindrical surfaces § and S, (Fig. 6) with generators parallel, respectively, 
to the axes OZ and OY, and directrices the flat curves Ty = f(x) (in the plane x = O)Jand Iy[z=—ky/p, f(x) 
(in the plane y = 0)}. 


Fig, 6, 


The line of intersection of Sy with the plane z = 0 is a straight line — the axis OY, while that of S, with § 
is a space curve I’; [y = f(x), 2=—Ky/ pp f(x)} 


pod | 

} 

| 
| 
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Let us examine the trihedral, open, infinite region Cy, bounded by the surfaces S, Sy, the plane z= 0 and the 


0], ry = f{(z),z= —=4(2) when x< 0], and the negative half of the axis 


oY (y < 0). Figure 6 shows part of the region Cy {O4, OG, OB;} . 


In the region Cy the derivatives x, y, Z do not change sign, If the trajectory X(t) = X[t, X(0)], at time t = 0, 
shares a common point X(0) with a boundary of Cy (the origin is excluded from the investigation), i.e., X(0) GC, /cx, 
X(0) #= {0, 0, O}, then X(t) leaves Cy strictly when t > 0.** 


Indeed, if X(0) = { x(0), y(9), 2(0)} lies on some edge of region Cy, for example on boundary S;, then the 
projection of the velocity vector 


X (0) ={2(0), y(0), 2(0)} on grad S, [x(0)], 0, —p,)* 


will be equal to 


X (0) grad S, = — [x (0)} {y (0) — f [x (0)]}} > 0; 


if X(0) lies on an edge of region Cy, for example on rf, then X(0) x grad Sy = 0, but X(0) gradS=z> 0 
(grad S = {-f [x(0)], 1, 0} ). 


Consequently, trajectory X(t), in both cases, leaves Cy whent> 0, The proof is analogous in other cases. 
From this it follows that closed trajectories cannot be contained in region C. 


Analogous examinations of an, in a sense, symmetrical (relative to the plane { OY, OZ} ) with t to 
C,;**** trihedral, open, infinite region Cy,, bounded by the surfaces S, Sy, the plane Z = 0 and edges >), rf”) 
and the positive half of the axis OY, permit one to conclude that for Cy; too, trajectories having, at t = 0, common 
points with its boundary, leave Cy; when t= 0. Moreover, region Cy also does not contain any closed trajectories, 


For what follows it is important to note that not a single trajectory X(t) = X[t, X(0)]J, X(0) F Cy, Cyy (X(0) 
&0)***** can adjoin the origin, i.e., | X(t)| = (G2 + —» when t —>+, because the straight line 
(30) is contained, together with the bases (end faces) F and F, of cylinder Cy, in the regions Cy and Cy. 


5. Proof of Existence of Cycle 


Let us examine, in the three-dimensional phase space E, the sum of the open regions Cy, Cyy, Cp (open 
cylinder, part 3); Q,=€C,UC,, UC», an additional, to this sum, closed region Q. and the intersection of 
the latter with a closed parallelipipedP(P + boundary, part 2) — a closed toroidal-shaped region w= QQ) P 
(Fig. 7). 

From Inequality (24) and the results of part 4 it follows that 

X(t) if X(0)€a, (31) 


where the open region w = w is a boundary. 


Moreover, let us prove that 
X (t) €@ whent > 2T, (32) 


if X(0)€o. 
For this it is sufficient, without losing generality, to justify (32) only for the points X(0) lying on the edge 
MNVW, where x(0) =£€, /(® < y (0) < r, \z (0) | “ r (Fig. 7), or on the edge MyNyV1W3, where 


* The slanted line / indicates the difference between two regions, C the closed, and C, the open regions, 
**I.e., not one trajectory can touch the boundaries of Cy. 

*** At the points of discontinuity f (x) = f (x + 0) or f ‘(x —0) (for simplicity we will assume f '(x) to be piece- 
wise continuous with ist kind discontinuities). 

****Since in the functional class A xf (x) > 0 when x =# 0. 

°*°**X(0) does not belong to the regions Cy and Cy. 
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| Y Trajectory X[t, X(0)], emerging from the point X(0) 


of the boundary MNVW leaves for suffi- 
ciently small value of t > 0 and, consequently, cannot 
return to this boundary for sufficiently large values of t. 


+) 

For a sequence ofinterval with respect to @, points 
{X,(O} » (ices, X_(0)€ w), for which X,(0)—> A(0) 
when n-*> o, it follows from the continuity of f(x) that for 
x each fixed number 6 = 0 


X[T 46, X,(0)] X(T + 6, X (0)] = X (0). 


From (31) it follows that, for each n, the point X[T +9 
+ Xp(0)] Consequently, X(0) € and since X(0) can- 
not belong to the boundaries MNVW and M,N,V;W,, then, 
of taking (24) into account, Relationship (31) is applicable 
Fig. 8. to this point. From this follows (32), which was what was 
required to prove, 


ies, 


Each trajectory X(t) = X[t, X(0)], for which X(0)€ w 
and X(0) does not lie in the plane z = 0, intersects the 
latter whent> 0, In the opposite case a contradiction 
arises because when y = z does not change sign, y is bounded 
(| y| < R, part 2 and, therefore, there exists lim y(t) = yoy 
to 
which determines the other boundaries, following from (12), 


(31) 


(32) 


In accordance with (32), the singular point { x,,, Yas 
Zept © W, which cannot be because W does not contain 
singular points, 


Let us examine the intersection of the sphere S (part 
2) with the plane z = 0 (Fig. 8). 


In the circle Sp, with radium R, there is a rectangle 
Ke My q (the intersection of the plane z = 0 with the paraliel- 
ipiped P) and the closed regions f9(B* H® PQ) and 


Fig. 9. 


> 2 , 
M 
Y 
ge 
LEA: 
8, 
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Hf PyQy). Inside these last two are the closed regions (Bg HPQ) and Gog (Byt (intersection of the 
toroidal-shaped, closed region | with the plane z = 0), 


If X(0) € weg, then, in accordance with part 2, the trajectory X(t) = X[t, X(0)] is Lagrange-stable in the 
positive sense, i.e., X(t)G@ S whent > 0, 


Moreover, for sufficiently small t = ty > 0, the representative point X(t) falls into that part of @, where z> 0 
and, consequently, the trajectory X(t) must intersect the plane z = 0 when t = tg> ty. This point of intersection 
Xq = X(tg) must lie in the region Q,, since into it enter all the trajectories for which z> 0 (part 4). 


Analogously, it can be proved that X(t), in leaving &, whent> ty, crosses Qy whent = T > ty, L.e., 
X(T) € Qos. 


The time T for one revolution of the representative point of the trajectory X(t) = X[t, X(0)], [X (0) € Dog] 
about the coordinate axis OX, for which X(T) € My, satisfies the inequality*® 


2n, 
(33 


where pg is the radius of the base of the cylinder Cy (part 3) while 


1 
Vnax = max |X (t)| = max (2? 4 y? + (R +f)? + + kf , 
= max {|/(®|, 7 


If one examines the section of sphere § made by the plane S,, passing thru the axis of the cylinder Cy, [the 
axis lies on the straight line (30] and the coordinate axis OY, then one observes sections analogous to those exa- 
mined earlier (Fig. 9). 


The region (B* * PyQ;') is analogous to region Mp3, while the closed convex** region (Bye 
is analogous to Wpq. 
There are no trajectories touching the region Q; because the scalar product of the gradient of plane 5 
grad S, = {4(D —), 0, 1} 
with the velocity vector of the representative point V = {, y, 2} 
V grad S, = }(D—) [y — (2)] — > 0 


for all points X(0) € Tg, 
dower ig he there is a corresponding point X(T) of the first intersection, when t = T > 0, of 
the trajectory X(t) = X[t, X(0)] with the region Mg, which can be written in short form as X(T) = F [X(0)}. 


In this expression, dai chow dee OR during the time > rmin» Where 
T min ! determined by (33). 


The established correspondence F is continuous and of a one-to-one nature, since the function f (x) satisfies 
the Lipshits conditions. Such correspondences are called topological or homomorphous, 


Region is made to correspond to region contained in ag. 


*For each t in the interval 0 < t < T, X(t) does not intersect Tips. 
* *1f Xq € then each point of the segment connecting them also belong to Xy + € 


where the positive quantities £y and €, satisfy £; + € = 1. This follows from the fact that f (x) is nonconvex 
while the section (P{Q, , Fig. 8) of the plane S, with the curved side of cylinder Cy is a straight segment. 
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In accordance with (32), for each point X(0)€ Wg , there will be required not more than N = 1 + E (2T/pin)® 
turns of the representative point of the trajectory X(t) =X[t, X(0)] about the OX axis, for the latter to arrive 
back, moreover into Wg This property can be written in short form as; FON) (wg) a 

The Bol'-Brouwer®* theorem can be applied to the convex figure we 

In the continuous mapping of a convex figure into itself, there exists a stationary point [22].°** The latter 
determines the cycle Cl. 


In this way, there exists a trajectory X(t) and the quantity T, > 0, called the period, which is such that « 
X(t + T,) = X(t) for allt. The question of whether there exist other cycles, besides Ci ,for System (12) remains 
open, 


6. Some Properties of the Cycle 


The cycle is called simple if the zeros of the component functions { x(t), y(t), 2(t)} are devisible by the zeros 
of their derivatives, of which there exists not more than two pairs in each period T,. 


Cycle Ci -simple, It is sufficient to show that the arc of Cl emerging from the point X(0) € wg , returns for 
the first time into w, to the point X(T.) = X(0) having made only one revolution about the OX axis, 4nd not N> 1, 
as it has already beer? established, The latter is guaranteed by the Levinson-Massera theorem [23, 18]; if in the 


topological (continuous and single-valued) mapping F of the flat figure ws, into itself PO) (as ) a for every 
n> N, then there is a stationary point. 


Cycle Ci-stable mode. The trajectory X[t, X(0)](= X(0) when t = 0) is called a'stable mode if, for each 
€> 0, there isa 5 =6(€)> 0 such that, for any pair of points X, and Xg, on the trajectory which, at time t = t, 
are separated by a distance |Xy—Xg| < 6, the distance | X(t, Xy) —X(t, Xg) | < € for any t(t ty, t= 


in the theory of dynamic systems, this “internal stability” property was first formulated by Franklin, who 


called it the S-property. The zero solution of Systern (12) X[t, 0) possesses this S-property and is therefore a stable 
mode which is called trivial, 


The cycle Cj X(t + T..)= X(t) is a nontrivial stable mode. In the opposite case there exists series of points 
on Chl nen: and {ty} n =13,..9<t, and the number €> 0 such that 


|X (tn, Xin) — X (tin: Xen) | >* 
and 


| Xin — Xon|—> O when n-> 00. . (34) 


The series contain converging sub-sequences and possessing 
the same property (34), with Jim Xin = lim Xon = Xo, limt, = to, <Te. 
Let us examine the limit 


L =lim|X (tn, Xin) — X Xen) |. 


On the one hand, in accordance with (34) L> 0, on the other hand, 


* The function (greater integer) Z (2T/T in) is equal to the smallest whole number enclosed in 2T/Tnin- 

**such a theorem was proven, in 1905, by the Latvian mathematician P, R, Bol’ and, independently of him, in 
1911 by the Dutch mathematician Brauwer [21]. 

*** The region @s_ was introduced into the investigation for the sake of clarity of discussion. It can be shown that 
the. region is homomorphous to a circle and, because “(@pq) Bol’Brauwer's theorem is applicable to 
it too, 

**** Periodic, nontrivial, stable modes also exist. 


) 
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L = lim [X Xin) — X (too, Xan)] + lim [X (th, Xin) — X (teor Xin)] + 


since the trajectory CI X[t, X(0)] is continuous with respect to X(0) for a fixed t, and is continuous with respect 
to t for a fixed X(0). The latter follows from the continuity of f(x). 


The obtained contradiction proves that (34) is impossible, and, therefore, CZ is a stable mode. The question 
of the "external" stability of the trajectory of the CZ cycle with respect to other transient modes remains open, 
as for example, the problem of the Liapunov stability of Cl , in particular its assumptotic stabilities and its regions 
of attraction [24]. 


The method set forth permits one to investigate the behavior of a three-dimensional, automatic control 
system, with one nonlinear control element, in phase space. However, the results can be extended to include 
three-dimensional and multi-dimensional systems with several different types of nonlinearity. 


The relationships presented permit one, in practical cases, to evaluate the adjustment limits of the transfer 
ratios of automatic machines, corresponding to different operating modes of the system, as well as the limiting 
values of the amplitude of possible oscillations. 


Particularly simple results are obtained for the piecewise-linear characteristics of a control element. 


The author thanks B. N. Petrov, V. V. Nemytskii, E. M. Vaisbord and Iu. P, Portnov-Sokolov for their valua- 
ble comments, which were taken into account in this work. 
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CONCERNING THE NOISE STABILITY OF 
PULSE-FREQUENCY TELEMETRY 


N. V. Pozin 
Abstract 


The real noise stability of pulse-frequency telemetering, 
which takes into account the basic methods of realizing means 
of transmitting information, is analyzed. The mean error and 
the mean-square error are taken as criteria of noise stability. 
The method of signal “discretization” is used in the analysis, 
Design formulae for the evaluation of noise stability are deduced, 


The noise stability of telemetering can be characterized quite completely by means of the mean and mean- 
square errors, In pulse-frequency systems utilizing pulse frequency meters,* these errors can be conveniently deter- 
mined by means of a mean number of false pulses, produced by interference in unit time.** This number character- 
izes the mean deviation of readings on the receiving apparatus from the actual value in the absence of noise, i.e., 
the mean error of telemetering. The mean-square error is determined from the mean square of deviations from the 
mean, 


Let us consider the problem of finding these errors by examining the low-frequency channel of a receiving 
device under the following conditions: 


1) the duration of the received pulses is equal to that of the spaces (i.e., the mark-space ratio of the pulses 
is two), with the pulse repetition frequency being between f, and fy; 


2) fluctuating (smooth) noise acts in the channel; 


3) the low-frequency channel of the receiving unit contains a low-pass filter with an upper cut-off frequency 
feo Foo = fg), @ bilateral clipper which cuts out from the incoming pulses a narrow band at a Uy (Fig. 1), a pulse 
frequency meter which responds to arbitrarily shaped pulses with an upper repetition frequency fy = fo. 


1, Analysis of Noise Stability by the Method of Signal "*Discretization" 


To analyze the interaction of a signal with noise let us divide the pulse repetition period 1/f (Fig. 1) into 
intervals of time At = 1/2f oo and let us call At the duration of a unit (elementary) time segrent, Let us desig- 
nate the number of unit time segments in a pulse or a space by 


*For example, capacitor frequency meters, etc. 

** We are concerned here with the dete>mination of errors, dependent on relatively strong noise, in the presence 
of which the initial distortion of the signal occurs at the expense of the division of pulses and spaces, Correlation 
methods of reception are not examined. The effect of weak, fluctuating noise, causing fluctuation of the leading 
edges, and introducing some very small mzan-square error component, is not taken into account. 
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where f is the pulse repetition frequency. 
Let us apply the term unit (elementary) distortion to the false space of duration At, formed at the output of 


the limiter, under the action of noise, in any part of the transmitted pulse, or to the false pulse, of the same dura- 
tion, in any part of a space, 


For simplicity, let us assume that the probability of unit distortion in the space P, and the pulse Pp are equal, 
i.e., 


Py = Pp =P. (2) 


Let us suppose that the probability Pof unit distortion is known.® Because the above defined quantity At 
represents the correlation interval of smooth noise at the output of the filter, the formation of unit distortions can 
be considered to be independent [1]. The probability P(m) 
of m out of n unit elements of time being distorted can be 


Y, es found by means of, for example, the binomial law of probs 
Pm = C™P™ (4 — (8) 
2nat where 
Fig. 1. Break-down of pulse repetition period into nl 
2n segments of time, each of duration At = 1/2fco. Cl“ 


is the number of combinations of n taken m at a time, 


Let us examine a half-period, If there are m unit distortions, then there can be different combinations of 
these, resulting, for example, only in a change in the shape of the transmitted pulse, if the end regions of the pulse 
are distorted (spaces), or in the break-down of the transmitted pulse into 2, 3, . . ., (m + 1) parts, each of which 
has a duration of not less than At. This break-down of the transmitted pulse is received by the receiver as the 
formation of 1, 2,..., m false pulses, 


To find separately the probabilities of the formation of various numbers of false pulses for m unit distortions, 
let us represent the probability P,, as a sum of probabilities P,,(k), where k is the number of false pulses;* * 


Pm = Pm (). (4) 
k=o 
Determining P,,(k) as 
Pm (k) = Am (k) P™ (1 =P)", (5) 
we represent the number of combinations Cf” as a sum of coefficients A,,(k): 


cm = An (6) 


*For the determination of P,, and P, (or P if P., = Ps) for some cases see the Appendix. 


**The method for calculating the mean mumber Kk of false pulses in a half-period, presented below, does not take 
into account the possible distortions in the edges of the pulse and space, due to which, distortions in the shape of 
the pulse (space) can become false pulses. However, taking into account the relatively small significance of the 
combinations not taken into account, this method is preferable from a point of view of simplicity and objectivity. 


From an examination of the combinations of Cf it is simple to find that the coefficient A,,(0), which 
shows how many combinations lead to a change in the shape of the transmitted signal, and the coefficients Ap (1) 
and A,,,(2), which indicate how many combinations lead to the formation of one and two false pulses, are deter- 
mined from the formulae 


Am (0) = m + 1, 
Am (1) = +1) m] [n — 


Arm (2) = + 1) m(m—1)} — 1) (n — m—2)}. 


It can be shown that, in the general case, the number of combinations which form k false pulses for m unit 
distortions out of a possible n, is determined by the formula 


k+1 (m + 1)! (n — m — 1)! (Ta) 


which is valid for all positive k < m = n~-{k + 1), as well as for k =-4 and m = n. 


The total probability of forming k false pulses, with the probabilities of forming from m*=k tom =n=(k+ 1) 
unit distortions taken into account, is equal to 


m=n—(k+1) 
Pm(h)- (8) 
m=k 


The mean number of false pulses, produced by noise, in a half-period, is equal to 


km 
k= P(1) 2P(2)+... + (hq) —P(—1) = (®) 


k=—1 
where k,,, is the maximum possible number of false pulses: 


for even n>2, 


kn= for odd 


—1 for n=1anin =2. 


Making use of (8) we obtain 
km m=n—(k+1) 
k= Ye (10) 
k=-1 m=k 


Under Condition (2), for each pulse repetition period there is a mean number of false pulses equal to 2k. 
It is obvious that 2k represents the mean, relative error in the readings on the pulse~frequency receiver. 
The absolute, mean, telemetering error, due to noise, at a frequency f is equal to 


* Let us note that the probability of losing the transmitted pulse or of the appearance of a number of false pulses, 
equal to k = —1 is the probability of distorting all m = n unit pulses, i.e., p(—1) = p”. 
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2fk. 

The value of Ayy effectively depends on frequency, since f appears as a multiplier and K is also a function 
of frequency, In this way, at different points in the frequency range used for telemetering there will be a different 
mean error due to noise, In telemetry, one is usually interested in some single value of error, reduced to the 
utilized range Af =f, —fy. 


Analogously to the accepted method of determination used in telemetry, the reduced mean error can be 
characterized thru the ratio of the maximum value of the absolute mean error A,, (max) Vet the range Af , to 
this range, i.e., 


bay ay (12) 


The reduced mean error can be defined in yet another way. Let us call the reduced mean error 54, the 
ratio of the value of the absolute mean error A,,, averaged over the range Af , to this range, i.e., 


bay = as 


Turning to the determination of the mean-square, telemetering error, let us examine the segment of time 
T, containing N periods of 1/f : Ne=Tf 


The time T is the persistence of the transmission.* 


The absolute mean-square error Ame, sq, 1s found from the square root of the dispersion of the distribution 
of false pulses over time T, The dispersion of this distribution is expressed thru the dispersion of the distribution 
over the time t = 1/2f, {.e., thru the mean value of the square of the deviation of the number of false pulses k 
from the mean number k. In the case of equal probability of the appearance of false pulses in any period, we have 


Arne 2N (k —k)? (15) 
The random quantity k can take on values~1, 0, 1, 2,.. ., ky, corresponding to probabilities P(—1), P(0), 
Let us expand (15) in the form: 


Kn 


Averaging (16) over the range Af, i.e., obtaining the value of 4... and relating it to the square of the 
maximum change in the parameter over time T, i.e., to T? Af*, we obtain the square of the reduced mean- 
square error 


A? 
(17) 


The above formulae for P(k) and k, obtained by the method of break-down of a signal into discrete amounts, 
are valid for integral values of n, If, whenn>2(f <f_./2), the points at which 4,, and A. 4q. are calcu- 
lated are rather closely spaced, then the region 1< n< 2(fog>f > fo 9/2), containing, in many cases, the whole 


*It is convenient to consider the time T to be the time for the setting-up of indications in the receiver, if the 
persistent sections in the transmitting system can be considered to be aperiodic. 
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operating frequency range, has only two boundary points of calculation, Clearly, the method of signal "discre- 
tization” does not yield a rigorous solution within this region.* However, one can make use of a graphical approx- 
imation and an experimental verification, Both these methods show that the relationships A,, = ¢(f) and ai 1" 
= » (f), im the region f./2to f are very nearly linear. Consequently, for the calculation of X, and 

sq, from Equations (11) and (16) for a series of values of f , with the subsequent determination of the arith- 
metic mean, as well as for the determination of the simplified equations derived below, one makes use of a linear 
approximation of the values of A,, and A? sq. between the calculation points f co/2 and f .,- 


2. Simplified Equations 


In practice, great value is attached to a simple expression for the mean number of false pulses ina half- 
period, representing the approximate equality 


k =(n-2) P, (18) 


obtained from the conditions k * A;(1)P and spread over all possible values of n, i.e., 00> n> 1, Substituting 
(18) into (11) and making use of (1), we obtain an important, in practice, expression for the mean absolute error 
due to noise: 


Day = & —2f) P. (19) 
“ay 
é 


Fig. 2. Distribution of the absolute mean error over the range 
of frequencies 0 to f .o. 


Equation (19) gives the value of the absolute error at any point on frequency range 0 tof... The distri- 
bution of the error over this range of frequencies is shown in Fig. 2. The physical explanation for the transition 
of A,, into the negative region when f > f,,, /2, lies in the fact that sufficiently short pulses can disappear 
completely in the presence of noise, but cannot be broken-down into parts, 


We integrate the modulus of (19) to obtain the mean over the range, assuming a uniform distribution of the 
telemetered parameter: 


1 


After integration and division by Af we obtain an expression for the mean, reduced error (13) in the form 


2 2 1 


*Such a solution can, in principle, be obtained by means of the distribution of the duration of the distortions, and 
the coefficient of correlation of the noise in a video-channel, 
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(19) 


, and 


In Expressions (21) and (26) the plus is used if the range Af includes the point f.9/2, 1.e., fy < f.,/2<fa, 
and the minus is used in the opposite case. 


If the utilized range Af is symmetrical about fog/2, if = AL, then from (2) we 
obtain 
Say =P. (22) 
Let us turn to the expression for the mean-square error, Assuming that the mean deviation k is small, 


taking into account the probability of the formation, in a half~ period, of only one false pulse (k = 1) when one 


unit pulse is distorted (m = 1), and setting P(k) = A,(1)P, we obtain an expression for the sum, which appears as 
a factor in (16), in the form 


— =|n—2|P. (23) 


Substituting (23) and (1) into (16) we find the value of the dispersion: 
Aine sq= T'| Lo — 2f| P- (24) 


Averaging over the range Af , we find 
2 
Aine ay (25) 
Integrating, and dividing the result by T*Af*, we obtain the square of the reduced mean-square error (17): 


If the utilized frequency range Af is symmetrical about f 2, then 
sq= TA}: (27) 


me Conclusions 


Besides the determination of errors, the obtained results aid in the evaluation and selection of a frequency 
range for pulse-frequency telemetering. In accordance with Equations (21) and (26), 5 gy and 6 me. sq, depend 
on the width of the range Af as well as on its position in the band 0 to fog of the low-pass filter. In Figs. 3 and 
4, the ratio of the middle frequency f,», in the range Af, to the cut-off frequency is plotted along the abscissae, 
while the values of the errors divided by multipliers, which are independent of Af, from the right sides of the cor- 
responding Expressions (21) and (26), are plotted along the ordinates, Curves are plotted for several values of 4 f, 
expressed in fractions of fg, and these characterize the movement of the range Af in the band 0 to foo. From 
Fig. 3 it can be seen that the mean error has a minimum value when the range Af is symmetrical about f /2, 
Figure 4 shows the same for the mean-square error. Moreover, from Fig. 4 we establish the fact that the mean- 
square error diminishes as Af approaches foo. 


To clarify the dependence of the noise-stability on the value of the range Af , for a constant relationship 
between Af and the filter's pass-band, it is necessary to analyze Equations (21) and (26), keeping in mind that P 
depends on the filter's pass-band (see Appendix). 
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Fig. 3. Variation of the mean error with displacement of the 
range Af in the frequency band 0 to f .,. 


Fig. 4, Variation of the mean-square error with displacement 
of the range Af in the frequency band 0 to f oo. 


From the analysis of Equation (26) for real values of noise (P< 0.1), it follows that, with increase in 4f , 
and a proportional increase in the band of frequencies in the channel, the ratio P/ Af increases, i.e., the numerator 
increases more rapidly, and therefore, the value of the mean-square error increases, An examination of Equation 
(22) yields an analogous result for the mean error. 


In this way, we can come to the following conclusions, 
1, It is expedient to locate the range Af as symmetrically as possible about f.,,/2. 


2. If the range Af and the practically arbitrary band of frequencies are connected thru some constant 
relationship, then, the narrower the range (or band) of frequencies, the greater is the noise~stability. This last 
condition opposes the condition for accuracy of telemetered transformations, i.e., the transformation of the 
measured parameter into a pulse repetition frequency and the reverse transformation at the receiving end. 


APPENDIX 


Comments Concerning the Initial Conditions of the Analysis 


1, It was assumed above that the probability P= P, = P, is known, For practical calculations it is necessary 
to indicate a method for finding P or Py and P,. To find these probabilities, it is necessary to know the laws 
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governing the noise distribution in the low-frequency channel (after the demodulator), 


a) Let us first examine the case when the fluctuating noise has a normal distribution at the input to the 
receiver, containing a low-pass filter, a bilateral clipper and a pulse frequency meter.* There is equal probabi- 


lity P for the formation of unit distortion in a pulse or in a space if the mark-space ratio of the pulses is equal to 
two and U, = U,/2. 


If the effective value of the noise voltage Uy in the filter pass-band, and the level of the bilateral clipping 
U, are known, then the probability P is found from the well-known formula 


P= P {Uy> = (28) 
where © is the symbol for the Laplace integral, 


z z* 
(2) = (29) 


The values of this integral are given in tables, 


b) Let us examine a receiver containing a band filter, a linear amplitude detector and the other circuits 
inherent in a pulse~frequency telemetering receiver — a bilateral clipper and a frequency meter, A high frequency 
signal, keyed amplitude-wise by telemetry pulses, arrives at the input of the receiver. 


If the noise has a normal distribution in the high frequency channel, then, after the linear detector, in the 
absence of a signal, the noise has a simple Rayleigh distribution, This is.a condition for finding P,, In the presence 
of signals, the joint distribution of a sinusoidal signal and noise is governed by a generalized Rayleigh distribution 
law, also called the Rice distribution, This is a condition for the determination of Py {1} 


In the case when the effective values of the noise voltage in the filter pass-band Uy, and the levels of the 


signal U,, and the bilateral clipping ®* are known, the probability P, of unit distortion in a space is found from 
the formula 


(oy (30) 


The probability P,, of unit distortion in the transmitted pulse is expressed by the equation 


2 
0 


where Iy is a zero-order Bessel function, 
Values of this integral are given in tables.* * * 


*In practice, this can be, for example, in the direct transmission of pulses without remodulation, Moreover, the 


distribution of fluctuation noise at the output of a linear detector is sometimes considered to be approximately 
normal for large signal-to-noise ratios, 

* *If the level of the bilateral clipping during the transmission of a space is Up, then, during the transmission of 
4 pulse, this level is calculated as U, —Up. 

***See, for example, the curve in 1] or [2] — which the probability P,, is found as a function of 


P 
for series of values of 


In the case considered here, one should use, in the equations in parts 1 and 2, the value of P equal to the 
arithmetic mean; 


(32) 


2. One of the initial conditions of the analysis is the existence, in the “video-pulse” channel of the receiver, 
a low-frequency filter® with cut-off frequency f.... The basic filter, limiting the frequency spectrum at the 
receiver input, is the sub-carrier, band filter with pass-band W. The interaction of noise and signal in this band 
appears as beats of different components of the signal with noise components, and the noise components between 
themselves, which form, after demodulation, a group of difference frequencies, These difference (low) frequencies 
pass together with the signal and distort the reception. 


The signai does not contain frequence components higher than W/2, and the main part of the noise spectrum 
appears in this band. Consequently, in the calculation of the noise stability of a telemetry receiver, having, at its 
input, a band filter with a pass-band W, one can assume that 


ko = (33) 
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A SINGLE-CYCLE MAGNETIC SHIFT REGISTER 


A. la. Artiukhin and V. Z. Khanin 


(Moscow) 


This paper includes an analysis of the functioning, approx- 
imate calculations and design and experimental! data on the 
simplest register circuit using one core per bit, The influence of 
variations in the length of the shifting pulses on register operation 


is investigated, 


Magnetic shift registers are widely used in the hardware of computers, automation devices, etc. Many 
register circuits employing magnetic materials with rectangular hysteresis loops are known [1-7]. These circuits 
provide not only for transmission of binary-coded signals, but also serve as a basis for constructing logical and 
controlling circuits, These circuits are simple, economical, reliable and have a long working life. 


Best known are the two-cycle registers, which require two cores and four diodes for each bit, With the pro- 
per choice of the turns ratio of the input and output windings in the connecting circuits, the number of diodes 
required can be reduced to two [9]. Two-cycle circuits are sufficiently stable in operation, Analyses of them 
can be found in works [2-4]. However, two-cycle circuits require two series of shifting pulses and a significant 
number of cores and diodes. Therefore, there is definite interest in single-cycle registers [7, 8], which require 


one core per bit and, correspondingly, one source of shifting pulses, 


There is given below an analysis of the functioning of a single-cycle register which explains the essential 
physical processes occurring during register operation, and approximate computations are also given. 


Principles of Register Operation 


Figure 1 shows the simplest circuit of a single-stage register. Each cell (stage) consists of a transformer 
with three windings (input winding Wy, shifting winding Wg and output winding Wg), a diode B, a condenser C, 
and resistor R. 
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Fig. 1. 
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The condenser in the connecting link temporarily remembers the state of the core in the form of its charge, 
Discharge occurs across the input winding after cessation of the shifting pulse, which allows the use of just one core 


per bit. 


It is necessary, for normal operation of the register, that when the n'th core is read, a signal be transmitted 
only to the following (n+ 1)’st core, In the circuit given, this condition is met. We consider the transmission of 
a "1" from the n'th core to the (n + 1)'st. Let the "1" be realized physically by the core's residual induction, B,, 
and let "0" be realized by -B.. The shifting pulse changes the core's magnetic polarity from state"1" to state 
"0". A current pulse flows through the connecting link, charging the condenser. The resistor R limits the cur- 
rent which is taken off by the input winding of the (n + 1)'st core, The incomplete blocking of this current is not 
important, since in the (n + 1)'st core, at the same time, the shifting pulse induces a field opposed to the field of 
the writing current, Charging of the capacitor continues until the n'th core has had its magnetic polarity comple- 
tely reversed, The shifting pulse must cease at this moment. The condenser begins to discharge across the input 
winding of the (n + 1)'st core, magnetizing it to state "1." With this there is induced, in the output winding of the 
(n + 1)'st core, an emf of such polarity that the diode prevents passage of current. The information stored in the 
n'th core is thus transmitted ahead only to the (n + 1)'st core, 


The reverse motion of the code is limited, since the condenser in the feedback path is shunted by the for- 
ward impedance of the diode, and can not be charged sufficiently. 


Analysis of Register Operation 


It is not possible to take into account all the factors which affect the operation of a register, Therefore, a 
number of common assumptions are made in our analysis: 1) the core material has a rectangular hysteresis loop; 
2) during the entire process of magnetic polarity reversal, the field acting on the core has aconstant strength, H,,;* 
3) the shifting current pulse has a rectangular shape; 4) the forward impedance of the diode during capacitor 
charging is constant, and equal to Rg; 5) the back impedance of the diode is small, and does not affect register 
operation; 6) the analysis does not take into account the induction, leakage or parasitic capacitance of the windings, 
or the capacitance of the diodes, 


In what follows, the following notation will be used: , is the residual magnetic flux in webers, H is the 
magnetic field strength in ampere-turns /meter, Ucg is the voltage on the capacitor at the moment when the 
shifting pulse terminates, t, is the charging time of the capacitor in seconds, J is the mean length of the core's 
magnetic lines of force in meters, try is the time taken to reverse the core's magnetic polarity when the capaci- 
tor discharges in seconds, S is the area of the core's lateral section in meters*, C is the capacity in farads, R is 
the resistance in the connecting link in ohms, Rg is the diode's forward resistance in ohms, tg is the total time for 
discharging the condenser in seconds and i, is the discharge current of the condenser in amperes. 


The transmission of information from one core to the next may be decomposed into two individual steps: 
discharge and charging of the condenser. 


Condenser discharge, The circuit for condenser discharge is given in Fig. 2. With a "1" read from the 
previous core, the condenser is charged up to the voltage Uc, at the moment the shifting pulse terminates and all 
the register’s cores have gone to state "0." The discharge current ig induces a change inthe core's magnetic 
field. By Kirchhoff's law, uc = Wyd@/dt + igR. 


The discharge current of the condenser is determined from the 


relationship: 
+ 
ry ie f Yopt ",) By our second assumption, the magnetizing field acting on the 
eS ou ie core, which we designate by H,,;, is constant throughout the period of 
Up magnetic polarity reversal, giving us 
Fig. 2. 


*H,, is the material's dynamic characteristic corresponding, in the static regimen, to the coercive force, Hy, is 
larger than the coercive force since it includes the effects of eddy currents and magnetic viscosity, Methods of 
measuring H,, as a function of the magnetization time are described in the literature (Cf., e.g., [4]). 


for 


all 


(1) 


The voltage on the condenser is 


uc = — (2) 

If we take into account that, by Kirchhoff's law, $(t) = , for t = 0 then, from (1) and (2), we get 
= ®,. (3) 
Figure 3 gives the graphs of (t) for various values of UC; The process of magnetic polarity reversal ter- 


minates at the moment t,,, when the magnetic flux in the core becomes equal to @,, 1.€., ®(tm) = ®,- 
By solving Equation (3) with respect to tm, we obtain 


— 


For complete reversal of the core's magnetic polarity, it is necessary that the following inequality hold: 


C — igR) 40,CW 


We find the least value of voltage on the condenser necessary for complete reversal of the core’s magnetic 
polarity from the formula: 


Ay, 
Uc, =2)/ + (4) 


With this, the time of magnetic polarity reversal equals 


_ 
ig H,,,! 


mi 


— RC. (5) 


If Uc '* Uc, the polarity reversal will be partial, 


If a "0" is read from the previous core then, at the moment when the shifting pulse terminates, a noise volt- 
age Uy will appear on the condenser. If, with this, the discharge current from the condenser induces in the 
input winding a field less than the coercive force, the flux in the succeeding core will not vary and, consequently, 
the noise will not accumulate in the register, With this, Uc, = iy Rs where ty, is the current corresponding to the 
coercive force, 


However, as shown by experiments, even for Uc, = gk there occurs no significant variation of the magnetic 
flux, This is related to the fact that the noise discharge current differs little from ic, but the time of its action 
on the core is limited, From the point of view of stability of register operation, it is desirable that the admissible 
value of Uc,/ Uc. be large. On the other hand, this would entail a large loss of the energy stored in the conden- 
ser when a “1” is transmitted (Uc, must not reverse the core’s magnetic polarity). In this case, the requirements 
of operational stability and economy are contradictory, It ts possible to take Uc fc, = 1/3, This magnitude is 
probably close to the optimum since, in this case, the energy loss will be about 10% of what is initially stored on 
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the condenser, and the noise voltage capable of being engen- 
@ : dered in the register may reach a value of one third the signal 


~6, 


voltage, which would not disrupt register operation. Taking 
Uc, /Uc, = 1/3, we obtain from (1), (2) and (4) that 
Any ®, 
Uc, = = R=3 . (6) 
Whence, 
R 
WwW, = H,,, ®, 
If the number of turns is chosen by this formula, then the 
Fig. 3. Te variation of flux as con- 
pe ab Pe sig of magnetic polarity reversal is found, from (5) and (6), to 
tm = 2RC. 


After magnetic polarity reversal, the input winding will 
have law impedance and discharge of the capacitor will occur 
exponentially, with time constant RC (Fig. 4). It is possible to 
consider that the condenser is completely discharged when the 
voltage on it reaches 0.1 Uc,, 1.e., (1/3)Uc, exp(tac/RC) = 


= 0.1Uc,. It follows that = 0.8RC, and the total discharge 
tn time equals tq = tr + tac = 2.8RC. 
td Se Condenser Charging, We now consider the most difficult 


; operating case, when a "1" is simultaneously read from all cores 
ft ~ pene ageing of the register. With this, in the input and output windings of the 
"6 n'th core there will be a flow of current opposing its magnetic 
polarity reversal and affected by the simultaneous polarity reversals of the preceding and succeeding cores (Fig. 5). 


The average value of the resulting field strength equals the algebraic sum of the field strengths acting on 
the core; 


It follows from Fig. 5 that 


is +i. 


We find, from the last relationships, that 


i cw, duo (H —H,,,)! 


(7) 


By Kirchhoff's law, 


i (Hz — H,,,) ! duc 
i+ WW, dad 
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Fig. 5. Condenser charging circuit, 
om = i,R, +U¢ am — => i,R — uc. 
By substituting in these equations the values for the currents ty and ig, we obtain the equation 


Cc R (W,/W;)*| (Ay _ Hq) l [R 


uc (1 + W,/Ws)* + 
the solution of which has the form: 
uc =E (1 — exp (8) 
where 
= — Ang) (R — [R+ (9) 
Ws (1 + W,/Ws) (1 + 


are defined by the boundary conditions that uc = 0 for t = 0 and duc/dt = 0 for t = a, 
The expressions for the currents ig and i, can be rewritten in the following form: 


CW,/W, t E (i + W,/Ws;) 
xp(—+)+ 


l = E+ Wi/Ws) (—<) 
R— ,/Ws i+WiW, P 


It is possible, knowing the laws by which the currents i, and iy vary in the connecting link and by which the 
voltage uc varies on the condenser, to determine the law defining the velocity with which the magnetic flux in 
the core varies during charging of the condenser: 


CRy(Wi/Ws)E E(Rg+ R) 


We now determine the time of condenser charging t, during which the core changes state, from ©, to —®y- 
At the instant t = t., @(t)=—@, and therefore, we obtain from the previous expression that 


E + t CRyWi/Ws | to 
20, = —E| Ws (~ 


The equation just obtained was not solved explicitly for te. But, from this equation and from (8), we obtain 
the ratio: 


R—RyWi/Ws CRqW;/Ws ] 
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tl LS) The right side of Equation (10) can be computed, By 
using the graph of (t./TY(1 —exp(t,/7)) = f(t,/T) (Fig. 6), it 
is possible to determine t,/r and, by computing T from For- 
/ mula (9), to determine the capacitor charging time te. 


28 

i é We return to Expression (8). For a constant magnitude of 
2 

22 


shifting pulse current, the voltage Uc, has a maximum value for 
VA some optimal ratio Wy/Ws. With this, the ampere-turns I,Wg, 
/ necessary for charging the condenser to Uc, will be a minimum, 


To find the optimal value of Wy/Ws analytically is quite 
difficult, An approximately optima! turns ratio may be found by 
determining the maximum of the function E =f (Ws) for a given 


16 Wy: 
uM Ws opt = Wy (1 + 2Rg /R) 
or 
06 12 16 20 26 26 (Wy/Ws) opt = + 2Rq). (11) 
Fig. 6. An exact value of (Wy/Ws) opt can be found by a series of 


tentative calculations, 


Reverse Signal Transmission, When a "1" is read from the n'‘th core, a voltage is induced in its input winding, 
As a result, there occurs a decrease of the condenser's charge (if a "1" was also read from the (n—1)'st core) or an 
overcharge of the condenser (if an isolated "1" were transmitted), In the latter case, after termination of the shift- 
ing pulse, there occurs a partial variation in the flux in the (n—1)'st core, which leads to an increase of noise at the 
following step. It has been shown experimentally that no further increase of noise occurs in the register, if the over- 
charge voltage, Ucy when an isolated "1" is read, is some five to ten times less than the voltage Uc, (Fig. 7). 


We now consider that R must be so that this last condition 
holds, 


The voltages Uc, @ and Uc, can be defined by means of 
Duhamel integrals: 


Fig. 7, For computing reverse signal trans- 
mission, 


Ue, (t) = hy (O) (t) +\ (A) 


own 


Uc, (t) = hs (0) us (t)+ 


where h4(t) = (1/(1 + R/Rg)) (1 -exp(-t/T )) and h(t); = -a/a + Rg/R)) (1 —exp (-t/T )) are transfer fonctions, and 
= CR/(1 + R/Rg). 


The voltages u;(t) and ug(t) are related by the expression, = ug(t) W/Ws. 
As the result, we get 


‘ =< — 


Ue, ‘ om 
i 7 \ uy (A) da. 
+h/ Ry 


By taking (11) into account, we find the desired formula: R = (a —2)Rg, where a is between 5 and 10. 


Influence of the Length of the Shifting Pulse on Register Operation 


We now consider the case, of practical interest when Ig > Igg and t, > t,, where Igy is the minimum shifting 
current Ig, for which, during time t,» ,the condenser is charged up to Uc,» and t,' is the time taken to reverse the 
core's magnetic polarity by current Ig. 


With this, there is complete read-out from the n'th core, At the moment of time th, Uc > Uc,- During the 
time At = t2 -t’,, the condenser partially discharges into the (n + 1)'st core, already shifted in aguas polarity, 
Two variations are possible; 


a) during the interval At, the condenser discharges to Uc < Uc,» and there is an attenuation of information 
in the register; 


b) during the interval At, the condenser discharges to Uc = Uc, (Fig. 8), with which normal operation of 
the register is assured. 


In order that uc(tg) = Uc,, it is necessary that 
At 
Uc = Uc, FG» (12) 


where U‘, is the voltage on the condenser at the instant t. 
We have, from Formula (8), 


It follows from Formula (9) that E = A(H, —H__), where 
i(R—R W,/W;,) 
A= H,=1,W,/l 


and H,, can be taken off the graph of H,, = f(t) for the instant t:. 
It follows from (12) and (13) that 


Fig. 8. Variation of condenser voltage 
in a forcing re n, 

g regime At 


l= (14) 


By assuming that At << tg, and using the notation 


Uc,! . 
we obtain 
I,= + yexp 


Thus, a variation of the shifting pulse length by At requires an increase in the current lg. During the inter- 
val At, the condenser ineffectively discharges, with time constant RC, through the writing winding of the (n + 1)'st 
core, which is saturated by the current I, during this time. 


To remove this flaw, more complex circuits are used, with LC delays in the connecting links (Fig. 9) or 
with subsidiary diodes, By (Fig. 10). 


ing, 
ft Uq (13) 
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From preceding 


stage 
To shifting 


pulse gene~ 
rator 
Fig. 9, Circuit for a one-cycle register with an LC lag in 
the connecting link. 
To follo 
From preced gan wing 
ing stage o shifting 
rator 
Blanking voltage 


Fig. 10, Circuit for a one-cycle register with a subsidiary diode By in 
the connecting link. 


In the latter case, a blanking voltage is applied during the time the shifting pulse is applied.® 


Register Design 


On the basis of the circuit analysis given above, it is possible to present a method for an approximate 
register design, 


Assumed as given are the material and dimensions of the cores, the diodes in the connecting links and also 
the speed of action. 


The sequence of design steps follows. 
1, The period between shifting pulses is T, = 1/f. 
2. Rand C are determined from the following formulae: 


a) R is between 5Rq and 10R,, b) ty + t, = Ty, where ty = 2.8RC. Tentatively assuming that t, = 0.5ty, we 
get that R < T;/4,2C. 


3, The time to reverse the core's magnetic polarity upon condenser discharge is determined from the for- 
mula try = 2RC. 


4, From the graph of H,, = f(t), Hy is determined for t = t,,. 
5. The residual magnetic flux is 4, = BrS. 


6, The voltage on the condenser necessary for complete reversal of the core's magnetic polarity is Uc, = 
3 Hy) 4. 


1. The number of turns on the input winding is Wy = Hy 
8. The number of turns on the output winding must satisfy the following inequality: Ws = Wy (1 + 2Ry/R). 
9. The time constant of condenser charging is T = C(R + Ry (W3/Ws))/(1 + W,/W;}. 

10, We compute the function: 


te R—RW,/Ws RyCW, /Ws 


*The circuit with the blanking voltage was first presented by V. A, Zhozhikashvili and K, G, Mitiushkin. 
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for- 


g/R). 


17, 
18, 


11, From the graph (Fig. 6, we determine a = t_/T, 
12. The time for charging the condenser is given by t, = Ta. 


13, We next determine the maximum speed of operation of the 
register, characterized by the quantity f, = 1/T, = 1/(2.8RC +t,). 
f_ must be greater than, or equal to, fy. Otherwise, it is necessary 
to decrease RC (Cf. step 2). 


Fig. 11. 


0 6 20% 26 32 sec 14 E= 


From the graph of H,, = f(t), Hyg is determined for t = t.. 


Hl = = 


E(i+W,/W;)W; 


The pulse power during one cycle is Pp = Uglg = 26, Wgl,/t.. 


The average power in one cycle is Pay = ptc/ Ty. 


Fig. 12. a) is the shifting pulse current I, and b) is 
the voltage on the condenser. 


Fig. 13, a) is the voltage on the condenser and b) is 
the current in the input winding. 
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Electrical Parameters of the Register 


Parameter Computed Data Experimental data 
Uc,» volts 6.5 5.2 
Wi 155 155 
Ws 185 185 
tc, microsecond 21.4 19 
We 200 200 
I,, milliamps 100 90 
tm, microsecond 48 45 
tq» microsecond 67 71 
Pp» watts 0.625 0.51 
Watts 0.134 0.1 
Speed in kilocycles 10 10° 


*The maximum speed equals 11 kilocycles, 


Fig. 15. 


In order to test the various assertions of the analysis and the computational method of design, an eight-digit 
register was designed and built. 


The core elements of the register took the form of annular rings, 9 mm in diameter and 3 mm high. The 
core material was alloy 50 NP of a 20 y strip thickness, The strip had four turns, The function H,, = f(t) is 
given on Fig. 11, The diodes used were germanium diodes, type D9G, with Rg = 20 ohms, The computed speed 
of operation was f = 10 kilocycles. 


For the generator of the shifting pulses, a 6P6 vacuum tube was used, The register worked in a ring circuit, 


It follows from the Table which gives the register's electrical parameters that the experimental and computed data 
coincide, 


Figures 12-15 show photographs of the voltage and current, illustrating the character of the basic physical 
processes during register operation, One photo (Fig. 13) confirms quite well one of the basic assumptions of the 
analysis and computations: the constancy of the condenser's discharge current during magnetic switching of the 
(n + 1)'st core. 
SUMMARY 


The analysis given of the simplest circuit for a one-cycle register permits its approximate design with 
an accuracy which is sufficient for practical purposes, 


The design method given can be generalized to the reading from two and more cores, 


It was shown that, in contradistinction to two-cycle registers, the one-cycle register is critically affected 
by variations in shifting pulse length. 


The reverse signal transmission occurring in the circuit, with its attendant shunting of the current to the 
(n + 1)’st core, when a "1" is read from the n‘th core, leads to a dependency on the structure of information in 
the register and of the form of the signal on the condenser. This might create difficulties in constructing logical 
circuits with single-cycle registers, 


The simplest circuit, considered here, can be recommended for use in ring counter circuits and in registers 
which are not too long. 


Received May 25, 1957 
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SCIENTIFIC SEMINAR ON PNEUMO-HYDRAULIC AUTOMATA 


From March 17, to March 19, 1959, at the Institute of Automation and Remote Control of the ANSSSR{Academy of 


Science of the USSR), there was held the Second All-Union Seminar-Meeting on Pneumo-Hydraulic Automata, 
organized in 1957, 


In attendance were 328 people from 63 organizations, 62 of the individuals representing organizations work- 
ing in the domains of pneumo- and hydro-automata in Leningrad, Baku, Kiev, Sverdlovsk, Khar'kov, Cheliabinsk, 
Dnepropetrovsk and many other cities, 


Taking part in the work of the seminar were specialists in questions of pneumo- and hydro-automata from 
other national democratic (Soviet block) countries: Prof. Lu luan'-tsin (KNR) (Chinese National Republic), 
Dzhozef Kveton and Ian Khampl (Czechoslovakian National Republic) and Britall and Maier (GDR)(German Demo- 
cratic Republic) and others. 


In opening the session, seminar leader Prof. M. A. Aizerman, presented the problems on which the seminar 
would work. 


First heard was the paper of Engineer V. A. Nikitin (Giprogaztopprom, Moscow), “The AUS pneumatic 
assembly system the basic automated complex in oil-refining.” The paper contained brief characterizations of the 
technological means for automatic inspection, regulation and control which are necessary for the automation of 
complex oil-refining processes, and also provided results of the experimental use of the AUS system. 


R. A. Auzan (TsNIIKA, Moscow), in her paper on the subject, “An investigation of the dynamic character- 
istics of the pneumatic regulators of the AUS system,” presented the design methodology and results of an experi- 
mental investigation, of the regulating blocks of the AUS system, 


Next, was read the paper of Engineer Ferner (GDR), “Pneumatic regulators and computing links.” The paper 
described analog pneumatic elements, developed in the GDR, which execute the computational operations which 
might be required for solving regulation problems, The system works at low pressures (0 to 100 mm of water), 
which only in extreme cases of necessity is increased to the range 0.2 to 1.0 atmospheres, Also considered in 
the paper was the superiority of a low-pressure system in comparison with those of the higher pressures which are 
usually employed, 


V. Britall (GDR) answered the questions raised by Ferner's paper. 


The evening session was opened by V. V. Volgin (Moscow Energetics Institute) with the paper, “Pneumo- 
hydraulic methods of implementing advancing.” The paper presented results of a theoretical investigation related 
to the use of pneumatic blocks for forward and reverse advancing in the regulation of technological processes by 
pneumatic regulators, The stability of single-loop automatic control systems containing advancing blocks was 
considered, and certain conclusions were given as to the influence of forward and reverse advancing on the 
quality of the transient response. 


L. L. Feigel'son (KB Tsvetmetavtomatika, Moscow) presented the communication, "A brief survey of the 
designs of pneumatic transducers in domestic and foreign production,” which was devoted to an analysis of the 
principal circuits and designs of pneumatic pressure transducers for the discharging and transmission of pressure. 
The designs of five domestic organizations and 23 foreign companies were considered. 


In his communication, entitled "The use of an extremum-regulator for the inspection and control of a 
reaction by certain chemical processes using the thermal effect," V. M. Dobkin (Moscow) related the work which 
he had carried out in conjunction with M. L. Kurskaia and Iu, O. Ostrovskii (IAT). In this work were considered 
variations of a scheme for regulating a process with two interacting components by means of an extremum 


regulator, and also a scheme for automatic thermotitration for determining the content of one of the products 
in the reactive mass, The results of experiments performed with laboratory equipment were given, 


E. M. Nadzhafov (IAT AN SSSR) in his paper on work done in collaboration with lu, I, Ivlichevol (IAT AN 
SSSR, Moscow) reported on the results of work in creating pneumatic universal multiplying-dividing devices, 
devices for extracting square roots and membrane~-less devices for adding continuously-varying pneumatic signals 
in the range from 0 to 1,0 atmospheres, 


A. A. Tal’ (IAT AN SSSR, Moscow), in the paper, "Pneumorelays and their applications in pneumatic discrete 


computing devices,” (work done jointly with T. K. Berends and N. V.Grishko) showed that discrete computing 
devices could be made on the basis of pneumatic apparatus, It was shown that the basic element of this apparatus 
is the universal pneumorelay, the circuit and miniaturized dimensions of which were developed in the IAT. The 
properties of this pneumorelay are such that it can be used for synthesizing all possible discrete pneumatic compu- 
ting devices. The circuits of four devices (a generator, a separator, a pulse counter and an interval indicator 

for a continuously-varying quantity), constructed of pneumorelays, were given, It was shown that it is possible to 
use pneumatic interval indicators in devices for the automatic load change-over of pneumatic regulators, 


L. A. Zalmanzon (IAT AN SSSR, Moscow), in the paper, “Aerodynamic oscillation generator (sine-wave 
generator)," reported on the development of a new type of oscillation generator with no moving parts, excitation 
of oscillation occurring only as a result of the interaction of air currents, Results of a theoretical and an experi- 
mental investigation of generator characteristics were reported, and the possibilities were cited of its use for 
obtaining the frequency characteristics of regulators and regulated objects, for driving various vibration equip- 
ments and in the solution of a host of other technical problems. 


L. A. Zalmanzon also reported on the development, done jointly with A. 1. Semikovoi (IAT AN SSSR, Mos- 
cow), of a multiplying-dividing device, constructed of nozzle~duct type elements, designed to function in 
pneumoautomata devices, With the introduction of an additional feedback path, this device can also implement 
the function of extracting square roots, 


V. N. Dmitriev (IAT AN SSSR, Moscow) reported on the piston follower pneumodrive, developed in the 
IAT AN SSSR, cited its basic parameters and gave its empirical static and dynamic characteristics, It was noted 
that in the drive distributed vents were used instead of valves. 


In the paper of E. V. Gerts (Institut Mashinovedemiia, AN SSSR, Moscow), “An investigation of resino- 
webbing membrane power drives," there were considered the existing formulae for determining the working stress 
or effective membrane area, and data were given from investigations made by the paper's author, A comparison 
of experimental and theoretical (computed) data was made . 


The morning session of March 18 began with the report by A. I. Semikovoi of work done jointly with L. A. 
Zalmanzon (IAT AN SSSR, Moscow) in the investigation of the characteris-ics of pneumatic chambers used as 


adders, The report gave data from an investigation into the influence of the compressibility of air on the accuracy 


of addition by means of pneumatic chambers, and ways were cited for increasing the accuracy of pneumatic 
chamber-adders, 


In the paper of the chief designer of the "Regula-vivoi? Dzhozef Kveton (Czechoslovakian National 
Republic), there was given a list of the pneumatic instruments released and recently developed by the factory, 
Two groups of pneumatic instruments have appeared: a basic group for the range of air pressure from 0.2 to 1.0 
atmospheres, and systems for controlling boilers, working in the pressure range from 0.7 to 3.5 atmospheres, 


The basic system was developed in the form of assemblies; right now, regulatory blocks with motion com~ 
pensation are being produced, and next year production will begin on regulatory blocks with power compensation, 
The speaker also reported on the executive mechanisms (membranes and pistons) produced in the factory, and on 
the work on the creation of new types of transducers. 


Dzhozef Kveton answered numerous questions from the conference participants, 


P. E, Baloban (VTI, Moscow), in his report on work done jointly with G. N, Makhan"kov, "Unified hydraulic 
regulating devices,” gave data on unified hydraulic regulators which wokkreliably when ordinary municipal water 
is used as the working agent. G, N. Makhan"kov (VTI, Moscow) provided additional information on experimental 
regulators which work on aqueduct waster and on turbine condensate, 
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V. Britall (GDR) gave a paper on the subject, "Hydraulic and composite hydraulic regulatory systems in the 
GDR." 


Two types of hydraulic regulators are produced in the GDR: regulators constructed on the assembly principle, 
for which systems of electrical remote task assignment have been developed, and normal hydraulic regulators with 
flow piping. Work is being done in combining hydraulic systems with low-pressure pneumatic systems. 


A number of measuring devices for hydraulic regulators have been developed (for measuring drops:under high 
static pressures, for working under the conditions of a viscous or agressive medium, etc.) 


A. M. Duel’ (Khar’kov), in his paper on “Certain questions as to the use of hydraulic regulatory apparatus in 
complex schemes of automated power plants,” considered the possibilities of using hydraulic auto-regulators for 
complex automation of the basic plants of thermal power stations. 


The evening session of March 18 opened with the paper, "Systems of coupled regulation,” of V. N. Veller 
(All-Union Heat Technology Institute, Moscow), The paper included comparisons of the various designs for 
coupled regulation of steam turbines. The disadvantages and advantages of each of the designs treated were 
mentioned, 


V. D. Mironov (VTI, Moscow), in the paper, “Electronic-hydraulic regulator, type EGP-T," spoke of the 
electronic-hydraulic regulator, based on an electronic regulation system and a unified series of hydraulic regulators, 
which was developed at VTI, All regulator elements have undergone lengthy tests under industrial conditions. 


Ian Khampl (Czechoslovakian National Republic), in his presentation, spoke of the work of the "Krzhyk- 
Smikhov" factory in the production of electro-hydraulic regulators. 


E. F, Alekseev, in the paper "On certain questions of the dynamics of rotating-piston drives as a result of 
experimental investigation," considered questions of the dynamics of electro-motor hydro-drives in relation to 
the characteristics of the former. 


According to the speaker, it has been experimentally established that oscillations in pressure and rotational 
speed of the output shaft, which arise under certain conditions, depend on particular design features on the machines 
and on the magnitude of the load on the drive. 


In the paper of V. A. Khokhlov (IAT AN SSSR, Moscow), “On one method of calculating the inertial load of 
hydraulic executive mechanisms with trottle control, by an analysis of servo system dynamics,” a method was 
presented for analyzing the quality of the transient responses in hydraulic and electro-hydraulic throttle-control 
servo systems, It wasshown that it is convenient, in engineering practice, to represent by an equivalent lag the 
inertial load connected with the input shaft of the hydro-motor in the analysis of servo system dynamics. A 
formula was provided for determining the magnitude of this lag. 


V. I. Gusakov gave a paper on the subject, "Increasing the speed of action of a throttle hydro-drive by 
stabilizing the pressure at the throttle aperture.” The speaker presented various methods of stabilizing the pres- 
sure in hydrosystems, 


In the report of V. P. Temnogo (IAT AN SSSR, Moscow), “Methods of increasing the gain of industrial 
hydraulic servos,” claims were made as to the efficacity of using hydraulic servo drives, both for the control of 
unitary processes and for complex automation, and an analysis of the basic parameters determining drive gain 
was provided, Six schemes of hydraulic servo drives were compared, 


V. M. Dvoretskii (IAT AN SSSR, Moscow), in his paper, “Hydraulic advancing blocks for general industrial 
use,” spoke of the principles of action and the characteristics of hydraulic advancing blocks of the compensated 
type. Data on their experimental investigation were given, 


In the opening session on March 19 B. L. Korobochkin (the Ordzhonikidze Stankozavod, Moscow), in his 
paper, “The automatic control of hydraulic transformer stands," spoke on work being done in creating new types 
of drives for the chief movement of stands ~ hydrodynamic transformers. The speaker stated that automatic 
control of the number of revolutions of a hydrodynamic transformer guarantees that rigidity of the latter will be 
obtained, which could significantly widen the domain of applicability of hydrodynamic transmission for driving 
machines and stands, providing a wide range of regulation with approximately constant power. 


ulators, 
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I. Z. Zaichenko (ENIMS, Moscow), in the paper, "On the dynamic stability of pneumatic and pneumo- 
hydraulic transmissions,” spoke of the results of investigations into the dynamic stability of volumetric hydrotrans- 
missions, Comparative data were given for the dynamic characteristics of hydrotransmissions, pneumotransmissions 
and combined hydropneumatic transmissions, 


V. N. Veller (Moscow) presented a paper, “An investigation of various types of piston hydraulic devices,” 
The paper gave several ways of decreasing the adjusting stresses on valves and pistons, An original design of a 
piston valve was discussed, | 


L. S. Bron (Moscow) made a report on the subject, “Hydraulic drives for power heads." Results of model 
research work on the creation of hydraulic systems of self-acting hydraulic power heads for stand assemblies were 
given. 


The paper of V. 1. Sherbakov (ENIMS, Moscow), "A new pneumoappliance for driving metal-working stands,” 
was devoted to designs for controlling pneumatic power drives for which automation of the machine's working 
cycle is guaranteed by just one pneumatic device, without use of electro-automatic elements, Also considered were 
the building of pneumatic terminal decouplers and pneumovalves, 


A. F, Arkhangel'skii, in his paper, "Ways to increase the power of hydraulic universal speed regulators within 
the previous dimensions,” spoke of experiments in increasing the power of oblique-plate hydrostatic transmissions, 
attained without any increase in physical dimensions, andrelated the results provided by these experiments, 


A. S. Shashkin (Avtomekhanicheskii Institut, Moscow), in reporting on the subject, "Nominal and structural 
nomenclature of hydrosystem elements,” mentioned the difficulties attendant on the study of hydromechanical 
semi-design schemes, and voiced his hope that standards of nomenclature for hydromechanisms and hydroapparatus 
would be developed and promulgated in the USSR. 


The evening session was given over to a discussion of the papers that had been heard, Taking part in the dis- 
cussion were E, M. Nadzhafov (IAT, Moscow), V. S, Prusenko (Moscow), P. M. Shanturin (Central Laboratory for 
Automation, Moscow), Iu. I, Ivlichev (IAT, Moscow), V. Britall (GDR), Io. 1. Nikolaev, S, A. Babushkin, V. A. Khokhlov 
(IAT, Moscow), L. S. Bron (Moscow), B. F. Stupak (The Leningrad Institute of “Sydostroitel’noi Promyshiennosti), 

Dzh. Kveton (Czechoslovakian National Republic) I, N. Kichin (IAT, Moscow, Lu, luan'-tsin (Chinese National 
Republic), S. M. Zasedatelev (MVTU, Moscow), S. A. Iushchenko (The "Manometr™ Factory, Moscow) and lan 
Khamp1 (Czechoslovakian National Republic). 


M. A. Aizerman summed up the work of the Second Meeting of the All-Union Seminar on Pneumo-Hydraulic 
Automata. He noted that the Proceedings of the First Meeting would be out in the near future, and spoke of the pre- 
parations for publishing the proceedings of the present meeting. 


A. 1, Semikova 


